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Preface 


As  the  former  leader  of  the  National  Aerospace  Plane 
(NASP)  Inlet  Team,  I  saw  that  boundary  layer  simulation  for 
inlet  testing  was  a  hurdle  that  needed  to  be  overcome  before 
the  inlet  flow  field  could  be  modeled  accurately.  Mass 
addition  seems  to  be  the  most  promising  method  to 
artificially  generate  a  boundary  layer.  At  the  time  of  this 
writing.  National  Aerospace  Plane  inlet  models  are  being 
tested  with  boundary  layer  injection  systems  which  were 
designed  empirically.  The  results  of  this  thesis  will  add  a 
theoretical  aspect  to  that  design  process. 

Software  discussed  in  this  report  may  be  obtained  from 
me  at  the  Ballistic  Systems  Division,  Norton  AFB,  CA.  The 
best  phone  number  that  I  can  offer  at  this  time  is  the  Norton 
AFB  operator's  assistance,  AV876-1110.  Software  may  also  be 
obtained  from  Mr  Don  Stava,  current  NASP  Inlet  Team  leader, 
WRDC/FIMM,  WPAFB  OH,  (513)  255-8508. 

I  would  like  to  express  my  appreciation  to  my  advisor, 
Lt.  Col.  Paul  King.  His  expertise  and  thoroughness  made  this 
effort  successful.  I  am  also  grateful  to  my  sponsor,  Don 
Stava,  for  his  time  and  experience. 
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Abstract 


A  computational  model  for  thickening  boundary  layers 
with  mass  addition  is  developed.  The  phenomena  of  uniform 
injection  into  a  two-dimensional  supersonic  stream  and 
subsequent  boundary  layer  growth  downstream  is  discussed. 
Analysis  of  the  injection  region  provides  the  thickness  of 
the  boundary  layer  just  aft  of  injection.  An  injection 
region  velocity  profile  is  then  used  to  approximate  the 
boundary  layer  profile  just  aft  of  injection  and  is  input 
into  a  finite-difference  boundary  layer  code.  Downstream 
profiles  and  thicknesses  are  calculated  and  compared  to 
experimental  results. 

The  computational  model  developed  here  provides  a  tool 
for  the  design  of  a  boundary  layer  generation  system  for 
hypersonic  engine  inlet  testing.  This  mass  addition  system 
is  needed  to  simulate  the  boundary  layer  developed  on  the 
forebody  of  hypersonic  vehicles^  An  example  is  discussed 
in  which  the  natural  boundary  layer  thickness  is 
increased  17  times. 


Introduction 


V.. 


"'The  United  States  Air  Force  is  currently  pursuing  a 
national  program  to  develop  the  technology  necessary  to  build 
a  single  staged  aircraft  capable  of  flying  from  conventional 
runways  into  orbit.  This  vehicle,  the  National  Aerospace 
Plane  (NASP),  will  use  airbreathing  propulsion  to  accelerate 
through  hypersonic  speeds  to  orbital  velocities.  In  order  to 
achieve  the  aerodynamic  and  propulsive  efficiencies  necessary 
to  accomplish  this  goal,  the  vehicle  forebody  must  be  used  as 
a  precompression  surface  for  the  engine  flowpath.  This 
results  in  a  thick  boundary  layer  which  must  be  bled-off  or 
ingested  by  the  engine.  Current  configurations  under 
consideration  for  NASP  will  have  boundary  layer  thicknesses 
at  the  inlet  cowl  station^Jon  the  order  of  20-30  percent  of 
the  cowl  height  (SJLava,  1989). 

^  ~  i  r%  . "" 

In  order  to  perform  a  wind  tunnel  test  of  an  engine  or 
inlet  concept,  large  plates  must  be  incorporated  ahead  of  the 
inlet  in  order  to  simulate  the  boundary  layer  expected  on  the 
vehicle.  Since  the  physical  size  of  the  model  is  limited  by 
the  test  facility,  the  size  of  the  inlet  is  forced  to  be  very 
small  in  order  to  accommodate  the  length  of  the  required 
extension  plates.  The  small  size  of  the  inlet  severely 
limits  the  amount  of  instrumentation  and  geometric  features 
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that  can  be  represented. 


In  order  to  demonstrate  this  difficulty  and  provide  a 
useful  example  for  later  comparison,  consider  the  boundary 
layer  thickness  on  a  hypersonic  vehicle  with  a  100  ft 
forebody  at  a  Mach  14,  120,000  ft  flight  condition  (Figure 

1).  For  a  simple,  single  ramp  forebody  geometry  the  99 
percent  boundary  layer  thickness  may  be  estimated  by 
( Schl icting ,  1951:638) 
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99 


0.37  x 


fD.  *  1 


“1/5 


(1) 


to  be  11.8  inches.  For  a  typical  inlet  cowl  height  of  about 
3  ft,  the  boundary  layer  is  32  percent  of  this  height.  In 
order  to  test  this  configuration,  a  Mach  10  facility  can  be 
used  with  the  boundary  layer  extension  plates  aligned 
parallel  to  the  free  stream  in  order  to  represent  the 
conditions  on  the  forebody  surface  (Figure  2).  A  suitable 
facility  is  the  Naval  Surface  Warfare  Center  Mach  10, 
Nitrogen  facility  which  has  a  large  test  rhombus  allowing  a 
model  over  6  ft  in  length  (Naval  Surface  Warfare  Center, 
1988).  A  6  ft  boundary  layer  extension  plate  in  this 
facility  would  generate  a  boundary  layer  height  of  0.50 
inches.  If  this  is  to  be  32  percent  of  the  cowl  height,  the 
cowl  height  is  only  1.5  inches!  Imagine  the  difficulty  trying 
to  simulate  internal  inlet  and  engine  geometry  or  incorporate 
instrumentation  in  this  small  size.  The  handicap  of  the 
boundary  layer  extension  plates  can  be  alleviated  if  a  means 
exists  to  artificially  thicken  the  boundary  layer. 
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One  such  method  is  the  use  of  mass  addition  (blowing)  to 
create  a  thickened  boundary  layer.  This  concept  would  be 
comprised  of  a  porous  plate  or  a  large  number  of  small  holes 
through  which  air  is  uniformly  injected  from  a  plenum  below 
the  plate  (Figure  3).  The  porous  section  would  be  followed 
by  a  length  of  plate  to  allow  the  boundary  layer  profile  to 
settle. 

Both  experimental  and  theoretical  work  has  been  done  to 
study  the  effect  of  mass  addition.  Experimental  efforts 
conducted  in  the  1960 's  concentrated  on  visualization  of  the 
streamline  between  the  injected  and  free  stream  fluid  and 
measurement  of  the  injection  surface  pressure  (Bott,  1967; 
Fernandez,  1968;  and  Hartunian,  1966).  Theoretical  analysis 
of  the  injection  region  has  yielded  streamlines,  velocity 
profiles  and  pressures  within  the  injection  region  (Lees, 
1968  and  Wallace,  1967  and  1968).  None  of  this  work  was 
oriented  towards  thickening  boundary  layers  and  did  not 
provide  information  on  the  flow  field  aft  of  injection.  One 
test  program  (Rozycki,  1968)  calibrate  a  mass  injection 
system  used  for  thickening  boundary  layers  in  a  ramjet  inlet 
test.  This  program  provided  velocity  profiles  downstream  of 
injection,  however,  a  generalized  model  for  predicting 
downstream  boundary  layer  thicknesses  and  profiles  is  still 
needed.  Thus,  the  purpose  of  this  thesis  is  to  understand 
the  injection  phenomena  and  provide  a  generalized  model  which 
can  be  used  to  design  boundary  layer  generator  systems  for 
hypersonic  engine  inlet  testing. 
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Three  regimes  exist  as  a  function  of  injected  mass  flux 
for  the  two-dimensional  problem  of  uniform  injection  into  a 
supersonic  stream.  The  first  occurs  when  the  injection  term 
(pQv  u^)  in  the  integral  momentum  equation  (below)  is 
comparable  to  or  smaller  than  the  skin  friction  term,  r 

o 

(Fernandez  and  Zukoski,  1968:1). 

,  f  Y  'i  w  f  nY  1  ctu 

r  =  —  |  pu2du  -  u  —  [  pudu  +  p  «  a  -  p  u  Y  — ®  ( 2 ) 

dx  o  J  dx  o  J  dx 

This  type  of  flow  results  in  a  perturbation  of  the  typical 

non-blowing  velocity  profiles  with  little  boundary  layer 

thickening . 

Second,  when  the  injection  term  is  large  compared  with 
the  skin  friction  term  and  the  injectant  mass  flux  is  small 
compared  to  the  free  stream  mass  flux,  the  injection  flow 
field  consists  of  the  following  regions:  a)  an  inviscid 
injectant  layer,  b)  a  shear  layer  interface  with  the  free 
stream  called  the  contact  surface,  and  c)  an  inviscid 
external  flow  field  (Wallace  and  Kemp  I,  1967:2)  (See  Figure 
4).  The  contact  surface  deflects  the  free  stream  flow  as  if 
it  were  a  wedge. 

Thirdly,  when  the  injectant  mass  flux  is  significant 
compared  to  the  free  stream,  the  incoming  boundary  layer  will 
completely  separate  resulting  in  an  irregular  flow  field. 
The  problem  addressed  in  this  report  is  solely  the  second 
regime  since  this  is  the  .  only  one  that  is  capable  of 
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generating  a  thickened  boundary  layer  with  a  steady  flow 
field  geometry. 

In  order  to  quantify  these  three  regimes,  two  parameters 

v 

are  useful.  The  first,  rj  =  —  Re1/2,  appears  in  Kays  and 

CD 

Crawford's  similarity  solution  for  the  laminar  incompressible 
boundary  layer  (Kays  and  Crawford,  1980:81).  A  value  of  q 
below  0.619  provides  boundary  layer  solutions  resulting  from 
wall  shear  (flow  regime  1).  Values  of  q  greater  than  0.619 
blow  the  boundary  layer  off  the  wall.  Kays  and  Crawford's 
analysis  pertains  to  incompressible  flow.  This  study 
addresses  the  compressible  problem  however  it  is  informative 
to  compare  the  values  of  17  for  this  study  to  r>  =  0.619  . 

This  study  deals  with  rj  on  the  order  of  1.  -  250  which  would 
indicate  that  the  the  flow  field  is  well  beyond  regime  1. 

The  second  parameter  is  the  ratio  of  injectant  mass  flux 

P 

to  free  stream  mass  flux,  X  =  — — -  .  If  the  integral 

od  p  ll  * 

K00  0D 

momentum  equation  is  nondimens ionali zed  by  free  stream 
dynamic  pressure,  -pa2  ,  the  injection  term  becomes  2-X 

2  CD  0D  OD 

T 

and  the  wall  shear  term  becomes  — 2 -  =  Cr  ,  the  skin 

12  0 

-pa 

2  OD  00 

friction  coefficient.  Although  dependent  on  Reynolds  Number, 
an  average  value  of  CrQ  without  mass  injection  is  about 
3xl0'3  -  4x10  3  (Lees  and  Chapkis,  1968:3).  This  study  deals 
with  X  from  0.005-0.08.  Since  the  nondimensionalized 

CD 

injection  term,  2-X  ,  is  from  3  to  50  times  Cr  ,  the 

oo  0 

injection  flow  field  type  is  regime  2  and  wall  shear  can  be 
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neglected  in  the  analysis. 

The  upper  limit  on  blowing  rate  for  regime  2  occurs  when 
the  contact  surface  deflection  angle,  5,  becomes  large  enough 
to  separate  the  incoming  boundary  layer  (if  present). 
Fernandez  and  Zukoski  experimentally  found  this  angle  to  be 
14  degrees  (Fernandez  and  Zukoski,  1968:3).  In  their  test, 
blowing  was  conducted  on  the  tunnel  wall  so  an  appreciable 
natural  boundary  layer  was  present  ahead  of  the  injection 
region.  Without  this  incoming  boundary  layer,  one  would 
expect  to  achieve  greater  deflection  angles  before  inducing 
separation  because  of  the  greater  momentum  near  the  wall. 
For  the  purpose  of  this  study,  it  is  advantageous  to  begin 
injection  at  the  leading  edge  of  the  extension  plate  in  order 
to  minimize  the  overall  length.  Since  no  incoming  boundary 
layer  is  present,  angles  greater  than  14  degrees  may  be 
possible.  However,  since  this  is  the  only  available  data 
point  for  a  separation  limit,  the  upper  limit  on  X^  for 
regime  2  flow  will  be  conservativel y  taken  to  be  the  value  of 
X  which  results  in  8  =  14°  . 

CD 

Injection  Region 

In  the  injection  region,  Wallace  and  Kemp  theoretically 
determined  the  geometry  of  the  streamlines  and  found  that 
they  become  closer  near  the  contact  surface  (See  Figure  4) 
(Wallace  and  Kemp  I,  1967:34,35).  This  indicates  that  the 
fluid  accelerates  as  it  is  forced  into  a  smaller  flow  area  by 
the  fluid  injected  beneath  it.  These  authors  also  found  that 
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the  pressure  is  constant  in  the  y-direction  and  thus  can  be 
matched  to  the  pressure  at  the  contact  surface  produced  by 
the  external  flow  (Wallace  and  Kemp  I,  1967:34).  Fernandez 
and  Zukoski 's  experimental  results  also  confirmed  the  absence 
of  any  appreciable  y-pressure  gradient  (Fernandez  and 
Zukoski ,  1968:3). 

Injection  may  be  accomplished  with  discrete  holes  or  a 

porous  material.  Injectant  velocities  are  a  function  of  the 

pressure  in  the  injection  region  unless  the  injection  holes 

are  choked.  Very  close  to  the  injection  surface,  large 

spatial  fluctuations  in  velocity  were  observed  by  Fernandez 

and  Zukoski  with  both  pitot  and  hot-wire  probes.  The 

fluctuations  decayed  rapidly  with  distance  from  the  surface, 

and  at  a  distance  of  0.1  inches  were  within  ±5  percent  of  the 

overall  mean  value.  Most  importantly,  the  overall  mean 

injection  velocity  at  the  plate  calculated  directly  from  the 

hot-wire  measurements  agreed  within  ±5  percent  with  the 

values  obtained  by  taking  the  measured  total  mass  flow  to  the 

plate  and  dividing  by  the  injection  region  density  and  the 

measured  plate  surface  area.  (Fernandez  and  Zukoski,  1968:2) 

Once  the  injection  region  pressure  is  determined,  this  allows 

the  injection  velocity  to  be  readily  calculated  by 

m  RT  m 

U  =  - L_  .  - L  — L_  .  (3) 

P.1.  p.  2i 

where  m  is  determined  by  the  choice  of  X  . 

i  ® 
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Contact  Surface 

The  injection  fluid  and  the  external  fluid  meet  only  in 
a  thin,  viscous,  mixing  layer  at  the  dividing  streamline 
called  the  contact  surface.  Static  pressure  and  streamline 
inclination  are  continuous  across  the  contact  surface; 
velocity  and  temperature  are  discontinuous.  In  experiments, 
Hartunian  and  Spencer  and  J.F.  Bott  utilized  a  chemical 
reaction  between  the  injected  and  free  stream  fluid  to 
visualize  the  contact  surface  as  shown  in  Figure  5.  A  major 
result  was  that  the  separation  streamline  was  observed  to  be 
straight  for  all  blowing  rates.  An  exception  to  this  is 
that,  near  the  leading  edge  of  the  blowing,  convex  curvature 
could  be  observed  (Hartunian  and  Spencer,  1966:1306).  The 
curvature  produces  a  normal  shock  in  this  small  region  which 
induces  a  high  pressure  on  the  plate  that  asymptotes  quickly 
to  a  constant,  greatly  reduced  value  (Wallace  and  Kemp  II, 
1968:11).  The  small  curvature  at  the  leading  edge  of  the 
straight  contact  surface  encompasses,  approximately,  the 
first  0.6  inches  of  the  injection  length,  1  (Fernandez  and 
Zuckoski ,  1968:2).  The  injection  lengths  considered  in  this 
report  range  from  3  to  24  inches.  The  small  leading  edge 
curvature  may  influence  the  solution  for  the  shorter 
injection  length  cases  but  the  effect  is  most  likely 
negligible  for  the  longer  lengths.  Slight  contact  surface 
curvature  is  also  present  at  the  end  of  the  injection  region 
where  the  contact  surface  turns  parallel  to  the  wall.  This 
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is  shown  in  Figure  6.  In  order  to  simplify  the  analysis,  the 
contact  surface  will  be  assumed  to  be  straight. 

The  contact  surface  is  actually  a  turbulent  mixing  layer 
of  finite  thickness  which  must  match  the  velocities  and 
temperatures  of  the  outer  and  inner  regions  across  its  width. 
A  gross  measure  of  the  contact  surface  thickness  can  be 
inferred  from  the  photograph  of  Hartunian  and  Spencer's  flow 
visualization  experiments  (Figure  5).  The  slip  line  grows 
approximately  at  a  rate  of  0.03  inches  of  width  per  inch  of 
length  or  3  percent.  This  thickness,  however,  is  assumed  to 
be  zero  for  the  injection  region  analysis  but  is  included  in 
the  later  boundary  layer  development  analysis. 

External  Flow  Field 

The  inviscid,  external,  flow  field  is  deflected  by  the 
contact  surface  as  if  it  were  a  wedge  of  equal  deflection 
angle.  Oblique  shock  relationships  establish  the  flow 
conditions  behind  the  shock  and  therefore  set  the  pressure  at 
the  contact  surface  streamline  and  through  the  injection 
region  down  to  the  wall. 

Aft  of  the  injection  region,  the  contact  surface 
streamline  turns  and  become  a  slip  line  parallel  to  the 
plate.  The  external  flow  passes  through  an  expansion  (Figure 
6).  Slip  line  angles  not  parallel  to  the  plate  are  not 
steady  state  conditions  for  the  flow  aft  of  the  injection 
region.  Such  a  situation  would  cause  the  flow  area  below  the 
contact  surface  streamline  to  change,  resulting  in  pressure 
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differences  above  and  below  the  slip  line  which  would  restore 
it  to  parallel.  The  pressure  in  regions  d  and  e  across  the 
slip  line  must  be  equal  since  there  is  no  flow  across  the 
slip  line. 

Mass  Flow  Balance 

The  key  to  solving  the  general  problem  of  the  injection 
region  geometry  lies  in  determining  what  phenomena  governs 
the  contact  surface  deflection  angle,  5.  Consider  the  mental 
experiment  where  injectant  is  instantaneously  introduced  into 

a  supersonic  stream  at  a  constant  X  .  All  of  the  injected 

00 

mass  flow  must  pass  through  station  e  (See  Figure  6)  so  the 
contact  surface  deflection  angle  will  increase  transiently 
until  the  injected  mass  flow  can  be  accommodated  at  e.  The 
conditions  everywhere  in  the  flow  field  are  changing  as  a 
function  of  5  and,  most  importantly,  the  changing  conditions 
at  station  e,  along  with  the  flow  height,  I  ,  determine  m  . 
Thus,  the  steady  state  contact  surface  deflection  angle  is 
established  when  m  =  m  This  fact  allows  solutions  for 

•  i 

the  contact  surface  deflection  angle  to  be  determined  as  a 
function  of  injectant  mass  flux  parameter  X^. 

The  following  analysis  lends  itself  to  two  separate 
discussions,  first,  the  modeling  of  the  injection  region  and 
second,  the  calculation  of  boundary  layer  development 
downstream. 
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The  goal  of  the  injection  region  analysis  is  to 
calculate  the  contact  surface  deflection  angle,  6,  and  the 
flow  field  conditions  at  each  station  as  a  function  of 
injection  mass  flux  parameter,  X  .  Referring  to  Figure  6, 

CD 

flow  field  stations  and  regions  are  established  as  follows. 

Region  a  is  the  free  stream  condition.  Region  b  is  the 

inviscid  region  aft  of  the  shock  above  the  contact  surface. 

Region  d  is  the  inviscid  region  aft  of  the  expansion  above 

the  slip  line.  Station  i  represents  the  injection  plane  and 

is  defined  at  the  plane  0.1  in.  above  the  injection  plate 
m. 

where  U  =  i-  as  previously  discussed  (Equation  3). 

i  P.At 

Station  c  is  at  the  end  of  the  injection  length.  Station  e 
is  defined  to  be  one  1  height  downstream  of  station  c  and  is 
assumed  to  be  the  point  at  which  the  irregular  flow  angles 
leaving  station  c  have  become  parallel  to  the  wall. 

A  control  volume  is  established  for  the  injection  region 
from  stations  i  to  e  (Figure  6).  The  control  volume  follows 
the  contact  surface,  the  slip  line  up  to  station  e  and  along 
the  wall  of  the  injection  region  0.1  in.  above  the  surface. 
A  few  simplifying  assumptions  are  necessary  to  render  the 
problem  tractable.  1)  The  shock  and  contact  surface  are 
assumed  to  be  straight  and  attached  at  the  leading  edge  of 
the  injection  region.  2)  Shear  forces  are  neglected  on  the 
contact  surface.  (Wall  shear  stress  was  previously  shown  to 


be  negligible  for  flow  regime  2).  3)  The  flow  height  1  is 

e 

assumed  to  equal  the  geometric  height  1=1  tan  <5  (not 

C  A 

depicted  in  Figure  6).  The  contact  surface  actually  begins 
to  turn  just  before  the  end  of  the  injection  region,  as 
depicted  in  Figure  6,  so  la  cannot  be  exactly  known.  Later 

results,  however,  show  that  the  solution  is  not  sensitive  to 

1 

x  • 
e 

The  following  analysis  of  the  two-dimensional  flow  field 
is  one  dimensional,  i.e.,  flow  field  conditions  represent 
averaged  quantities.  For  the  control  volume,  the  continuity, 
axial  momentum  and  the  energy  equations  are  given  as 
Equations  4-6.  Each  mass  flow  term  is  per  unit  width, 
continuity: 

m  =  p  U  1  (4) 

i  0  e  • 

axial  momentum: 

p  (J  tan6)  -pl=(pUl)u-mU  cos6  (5) 

bi  aeaaaaii  i 

energy : 

m.(h. +uf/2)  =  (p  U  1  )(h  +U2/2)  (6) 

By  making  the  following  substitutions:  perfect  gas  law  for  p, 
1  =  1  =  1  -tan6,  and  h  =  C  T,  the  three  conservation 

a  c  l  p 

equations  become, 
continuity: 

P 

m  =  - 2-  U  I  (7) 

1  R  T  *  * 

e 

axial  momentum: 

P  , 

(p  -  p  )1  tan6  =  - —  U  1  tan5  -  m  0  cos5  (8) 

b  a  i  o  m  •  *  ‘  ‘ 
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energy : 


m  (C  T  +  U2/ 2)  =  - —  U  1  tanS  (C  T  +U2/2)  (9) 

i  p  i  i  _  _  e  i  p  q  q 

R  T 

e 

Setting  pe  =  pd  in  Eq.  9  and  solving  for  T  , 

e 

p  U3I  tanS 

T=  - : - V-1 -  (10) 

2Rm  (C  T  +  u  /2)  -  2p  U  C  1  tanS 

ipi  i  d  •  p  i 

Thus  T  =  T  (U  ,  5,  knowns).  Substituting  Eq.  10  into  Eq.  8 

a  e  e 

and  solving  for  U  , 

8 

2m  (C  T  +  U  2 / 2) 

U  =  - - - S - - - : -  (11) 

®  (p  -p  +  2p  C  /R)I  tanS  +  m  U  cosS 

b  d  dp  i  i  i  i 

Thus,  U^=  U^(6,  knowns)  only. 

The  test  facility  defines  the  free  stream  conditions, 
and  Tt>.  The  designer  specifies  the  injection  system 
parameters:  length  of  the  injection  region,  1  ,  injection 

i 

angle,  5  ,  and  the  value  of  X  . 

i  0D 

The  following  method  is  then  used  to  find  5.  This 
method  is  coded  in  the  Fortran  computer  program 
INJECT  (Appendix  B)  and  is  presented  in  flowchart  form  in 
Appendix  A.  The  shock  angle  8  is  first  arbitrarily  assigned 
a  value  of  10  degrees  greater  than  the  Mach  wave  angle.  A 
unique  8  can  then  be  calculated  (as  opposed  to  the  two 
solutions  for  0  as  a  function  of  5)  from  oblique  shock 
theory.  The  external  flow  field  conditions  (regions  b  and  d) 
are  then  calculated  using  oblique  shock  and  expansion  theory. 
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Since  dp/dy  a  0  within  the  injection  region,  p  is  set  equal 
to  p  .  U  is  then  calculated  from  Equation  3,  repeated  here, 

m  RT  m 


0  = 

i 


p.  i 


i  i 


P.  1  . 


The  pressure  at  station  e  (p  )  is  set  equal  to  p 


and  T  are  calculated  from  Equations  11  and  10, 

e 

The  mass  flow  at  station  e  is  then  calculated  by 

p  U  1 . tan6 


(3) 

and  then  U 

d  e 

respectively . 


m  =  p  U  2  « 

0  0  0  0 


RT 


(12) 


m  -m 

If  me  equals  mi  (within  a  tolerance  of  — - - —  <  10-3),  6  is 

m 

i 

found.  If  m  is  not  within  tolerance,  0  is  updated  using  a 
Newton-Raphson  iterative  method  until  5  is  found.  The  value 
m  -m 

of  — l— — S—  has  only  one  root  as  a  function  of  0  which 

m 

i 

facilitates  iteration. 


Resul ts 

Lambda  was  varied  between  values  of  0.005  and  0.075  for 
Mach  2.6  free  stream  conditions  and  a  3.4  in.  long  injection 
region.  These  limits  on  ensure  that  the  second  regime  of 
injection  flow  fields  is  obtained.  The  free  stream  and 
injection  conditions  used  are  the  same  as  the  tests  conducted 
by  Fernandez  and  Zukoski  to  allow  for  the  comparison  of 
results.  Results  of  the  analysis  are  summarized  in  Table  2 
and  graphically  presented  in  Figures  7  and  8. 

Figure  7  presents  contact  surface  deflection  angle  8  and 
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injection  pressure  and  velocity  as  a  function  of  X  .  The 

® 

contact  surface  deflection  angle  increases  steadily  with  x 

CD 

and  reaches  the  14  degree  separation  limit  at  a  value  of  X  = 

© 

0.07.  Pressure  in  the  injection  region  increases  with 

injected  mass  flow  and  ranges  from  1.2-2. 5  times  the  free 

stream  static  pressure.  The  injection  velocity  increases 

with  X^  and  ranges  from  only  1  to  7  percent  of  the  free 

stream  velocity.  This  indicates  that  the  injectant  momentum 

will  be  small  relative  to  the  free  stream  momentum. 

Figure  8  presents  the  conditions  at  station  e  as  a 

function  of  X  .  Since  6  is  relatively  small,  the  loss  in 
© 

total  pressure  due  to  the  shock  is  small.  Thus,  the  equal 
angle  compression  and  expansion  results  in  p  being  nearly 

d 

equal  to  the  free  stream  value,  as  can  be  seen  by  the 

straight  line  in  the  figure.  As  X  increases,  the  velocity 

© 

at  station  e  increases  to  accommodate  the  additional  mass 
flow  as  evidenced  by  the  increase  in  M  .  The  flow 
accelerates  from  the  injection  region  to  station  e.  Theratio 
U  /  U  ranges  from  24.3  to  7.75  for  X  =  0.005  to  0.075, 

•  i  © 

respectively.  T  decreases  with  increasing  X  because  U  is 

•  ©  • 

increasing  and  the  thermal  energy  has  been  converted  to 
kinetic  energy. 

The  experimental  results  of  Fernandez  and  Zukoski  appear 
to  be  the  only  source  of  empirical  data  for  5  as  a  function 
of  X  with  a  similar  geometry.  Their  Mach  2.6  tests  used  a 

CD 

3.4  in.  long  porous  plate  with  side  walls  projecting  into  the 
stream  to  preserve  two-dimensionality.  The  contact  surface 
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deflection  angle  was  determined  from  schlieren  shock  angles 
and  oblique  shock  theory.  Contact  surface  deflection  angle 
versus  X  for  the  test  data  is  shown  by  the  solid  line  on 

GD 

Figure  9.  Results  of  the  analysis  are  shown  by  the  dashed 
lines  with  the  baseline  case  of  1  =1  -tanS  as  the  center  of 

q  i 

the  three  dashed  curves.  The  theoretical  values  provide  good 
results  at  lower  values  of  X  .  At  higher  values  of  X  ,  the 

CD  00 

theoretical  analysis  predicts  lower  values  than  the 
experimental  results.  In  order  to  find  a  possible 

explanation  for  this  difference,  the  free  stream  and 
injection  conditions  were  varied  by  ±10  percent  to  determine 
the  sensitivity  of  the  solution  to  these  parameters.  If  the 
solution  was  sensitive  to  one  or  more  parameters,  an 
uncertainty  in  the  given  value  may  be  significant  enough  to 
cause  the  difference.  The  results  of  the  parametric 
variations  are  presented  in  Figures  9-14. 

Figure  9  shows  the  effect  of  a  ±10  percent  change  in  1  . 

e 

One  assumption  in  the  analysis  was  that  1  equaled  the 
geometric  height  1  -tan6.  This  variation  has  little  effect 

i 

on  the  solution,  especially  at  the  upper  values  of  X^  and 
thus  validates  the  assumption. 

The  ±10  percent  variation  of  the  free  stream  Mach 
number,  Figure  10,  has  a  moderate  impact  on  the  solution. 
The  free  stream  Mach  number  affects  the  oblique  shock 
pressure  ratios  and  the  value  of  m^  calculated  from  the  given 
value  of  X  .  Ten  percent  variation  in  free  stream 

CD 

temperature,  Figure  11,  has  a  limited  effect.  The  effect  of 
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varying  would  be  manifest  by  a  change  in  the  calculated 
value  of  m.  .  The  variation  of  free  stream  pressure,  Figure 
12,  has  no  effect  on  the  solution  which  demonstrates  that  the 
solution  is  a  function  of  pressure  ratios  and  not  absolute 
pressures.  For  a  given  X  and  M  ,  the  pressure  ratios 

CD  a 

throughout  the  flow  field  are  constant  and  independent  of 
variation  in  the  magnitude  of  p  . 

a 

Variations  in  injection  conditions  also  had  little 
effect.  Variation  of  T  directly  affects  the  calculated 

i 

value  of  m  but  did  not  significantly  impact  the  values  of  6 

i 

as  can  be  seen  in  Figure  13.  A  significant  result  is  that 
the  variation  of  the  injection  angle,  6  ,  had  almost  no 

i 

impact  (Figure  14).  The  ratio  of  injectant  momentum  flux 
to  free  stream  momentum  flux  ranges  from  10-4  to  5x10  3. 
Since  the  injected  momentum  is  small  compared  to  the 
free-stream,  the  direction  of  the  injected  momentum  has 
little  affect  on  the  flow  field.  This  indicates  that 
injected  momentum  plays  little  part  in  establishing  the 
contact  surface  deflection  angle  and  supports  the  basic 
premise  that  conservation  of  mass  drives  this  phenomena. 

In  summary,  the  solution  of  5  as  a  function  of  X^  is  not 
sensitive  to  the  free  stream  or  injection  conditions.  Thus, 
an  uncertainty  in  the  value  of  these  conditions  probably  does 
not  explain  the  difference  in  experimental  and  theoretical 
results  at  the  higher  values  of  X  . 

CD 

One  possible  explanation  of  this  effect  is  the  viscous 
interaction  at  the  contact  surface.  Viscous  forces  would 


create  a  thickened  mixing  layer  at  the  contact  surface  with 
lower  velocities  than  the  inviscid  region  b  velocity.  This 
would  create  a  displacement  thickness  that  would  effectively 
increase  6  and  the  shock  angle.  J.F.  Bott  conducted  similar 
blowing  experiments  and  found  that  the  effect  of  the 
interaction  is  to  push  the  shock  wave  to  a  larger  angle  than 
would  exist  without  the  interaction  (Bott,  1967:13).  Figure 
15  is  a  plot  of  shock  angle  as  a  function  of  effective  wedge 
angle  for  Bott's  experimental  results  and  for  oblique  shock 
theory.  In  all  cases,  the  experimental  shock  angles  were 
above  the  theoretical . 

A  second  possible  explanation  of  the  difference  between 
the  experimental  and  theoretical  results  is  the  curved 
contact  surface  at  the  leading  edge  of  the  injection  region. 
The  curvature  displaces  the  entire  contact  surface  in  the 
y-direction  before  the  linear  contact  surface  begins.  By 
neglecting  the  curvature,  this  displacement  is  not  accounted 
for  in  the  analysis. 

A  third  explanation  is  perhaps  that  some  anomaly  exists 
in  the  test.  The  side  walls,  for  example,  will  develop  small 
boundary  layers  that  will  displace  flow  within  the  injection 
region.  The  contact  surface  would  be  pushed  to  a  slightly 
higher  angle  to  accommodate  the  fluid  displaced  by  the  side 
wall  boundary  layers.  The  actual  effect  is  most  likely  a 
combination  of  the  three  possible  explanations.  The 
theoretical  model  however,  does  provide  a  reasonable  model  of 
the  injection  phenomena. 
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IV .  Boundary  Layer  Development 


Analysis 

For  the  region  aft  of  the  injection,  the  goal  is  to 
incorporate  the  results  of  the  injection  region  analysis  into 
a  boundary  layer  code  and  determine  boundary  layer  heights 
and  profiles  at  a  desired  location  downstream.  Station  e 
serves  as  the  starting  location  for  the  development  of  the 


boundary 

layer. 

The 

one-dimensional  analysis 

of 

the 

injection 

region 

yielded 

a  single  averaged  value 

for 

the 

velocity  and  temperature  at  station  e.  This  uniform  value 
was  input  into  a  finite-difference,  boundary  layer  code  as  a 
uniform  profile.  The  velocity  and  temperature  profiles  at 
station  e  are  highly  non-uniform  however  and  a  more  accurate 
downstream  profile  may  be  obtained  by  using  a  more  realistic 
representation  of  the  starting  profile.  The  actual  flow 
field  in  region  e  is  extremely  complex  and  cannot  be 
discerned  without  a  computational  analysis  using  the  full 
Navier-Stokes  equations.  However,  we  can  obtain  a  form  of 
the  velocity  profile  from  previous  work  describing  the  flow 
field  within  the  injection  region. 

Wallace  and  Kemp  mathematically  analyzed  the  flow  field 
within  the  injection  region  for  incompressible  flow  using  a 
stream  function  approach.  Their  result  for  the  velocity 
profile  is  (Wallace  and  Kemp,  1967:33) 

u  =  V  -(1-f  -2c  Inf ) 1/2  (13) 
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where  u  is  the  x-component  of  velocity,  V  is  the  injection 

W 

velocity  at  the  wall,  c  =  — and  f  is  a  transformation  of 

252 

y  given  by  y  =  I  (1-f).  The  value  of  f  is  defined  as  unity 

e 

at  the  wall  and  zero  at  the  contact  surface  and  is  taken  to 
vary  linearly  between  the  two.  A  plot  of  Equation  13  (with  f 
transformed  to  y)  is  shown  in  Figure  16.  At  the  wall  u  =  0 
and  at  the  contact  surface  u  approaches  infinity.  Although 
Equation  13  was  derived  for  incompressible  flow,  the  authors 
compared  the  incompressibl e  streamlines  to  those  for  the 
compressible  case  and  found  them  to  be  almost 

indistinguishable  (Wallace  and  Kemp,  1967:33). 

Since  the  contact  surface  slip  line  is  a  viscous  mixing 
region  of  finite  thickness,  the  actual  velocity  never 
approaches  infinity.  As  previously  discussed,  the  contact 
surface  thickness  at  any  given  x-location  is  approximately  3 
percent  of  the  length.  The  lower  half  of  this  thickness 
would  occupy  the  upper  portion  of  the  injection  flow  field 
region  governed  by  the  'f’  velocity  function  (Equation  13). 
The  lower  edge  of  the  contact  surface  thickness  defines  the 
upper  most  value  of  f  (f  )  (See  Figure  16).  For  a  given 

up 

value  of  V  and  f  ,  the  velocity  at  f  (u  )  is  known  from 

v  up  up  up 

Equation  13.  The  velocity  profile  across  the  contact  surface 

thickness  must  match  the  velocity  (u  )  and  the  slope  (du/df) 

up 

at  the  lower  edge  and  the  region  b  velocity  (U  )  with  slope 

b 

du/dy=0  at  the  upper  edge.  A  third  order  polynomial  is  used 
to  model  the  velocity  profile  through  the  contact  surface 
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thickness  and  these  four  conditions  are  used  to  determine  the 
coefficients  of  each  term. 

The  total  enthalpy  in  the  injection  region  is  given  by 
h  =  C  T  +  U2/ 2  and  is  constant  along  a  streamline.  The 

t  i  pi  i 

temperature  within  the  injection  region  is  determined  by  T(f) 
=  (h  -  u2(f )/2)/C  .  The  temperature  through  the  contact 

t  i  p 

surface  thickness  is  also  modeled  by  a  third-order  polynomial 
with  the  temperature  and  temperature  gradients  matched  at  the 
edges.  This  is  presented  graphically  in  Figure  17. 

The  'f'  function  and  third-order  polynomial  model 
applies  to  the  profiles  within  the  injection  region;  however, 
velocity  and  temperature  profiles  are  required  at  station  e 
in  order  to  begin  the  boundary  layer  analysis.  The  injected 
fluid  accelerates  from  station  c  to  e  but  it  would  seem 
reasonable  to  assume  that  the  form  of  the  velocity  profile  at 
station  e  is  similar  to  that  in  the  injection  region.  Thus, 
this  model  is  applied  at  station  e.  Since  station  e  is  aft 
of  the  injection  region,  no  longer  has  physical 

significance  but  is  only  a  parameter  which  allows  the  model 
to  be  applied.  In  order  to  solve  for  the  profile  at  station 
e,  V"  is  assumed  and  u(f)  and  T(f)  in  the  lower  region  and 
u(y)  and  T(y)  through  the  slip  line  are  determined.  The  mass 
flow  is  then  calculated  by  numerically  integrating  the 
profiles  as  follows: 
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m  =  m  +  m 

f  profile  slipline 


nfuP  p  1  u(f)  ie  p  u(y) 

=  r  __2_S -  df  +  I  —2 -  dy  (14) 

Ji  RT(f)  yup  RT(y) 

The  mass  flow  is  then  compared  to  m*  and  V  is  iterated  upon 
until  the  mass  flow  equals  m  (=m  ).  The  resulting  profile 
is  then  input  as  the  starting  profile  in  the  boundary  layer 
code . 

This  profile,  by  definition,  conserves  mass  at  station  e. 
However  in  order  to  be  valid,  momentum  and  energy  must  also 
be  representative  of  the  values  calculated  from  the  averaged, 
1-D  quantities.  Momentum  and  energy  are  numerically 
integrated  at  station  e  in  the  program  INJECT.  For  the  Mach 
2.6  case  previously  discussed  and  for  a  Mach  4.38  and  Mach  10 
case  to  be  discussed,  the  difference  between  the  1-D  and  ’f' 
profile  momentum  were  10.9%,  11.3%  and  15.4%  respectively. 
The  difference  between  the  1-D  and  'f  profile  energy  was 
0.0%,  0.3%,  and  12.6%  respectively.  Although  it  is  not 
possible  to  accurately  conserve  all  three  quantities,  these 
differences  are  low.  Thus,  the  'f'  profile  satisfactorily 
represents  the  1-D  momentum  and  energy  at  station  e. 

The  station  e  velocity  and  temperature  profiles  are  now 
input  into  a  boundary  layer  computer  code.  The  code  used  is 
TEXSTAN,  an  updated  version  of  STAN5  by  Kays  and  Crawford. 
The  code  provides  a  finite-difference  solution  of  the 
boundary  layer  equations  and  is  used  here  to  solve  the 
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problem  of  turbulent  boundary  layer  growth  on  a  flat  plate 
with  zero  pressure  gradient.  The  Reynolds  number  at  station  e 
was  calculated  for  the  Mach  2.6,  4.38  and  10.0  cases  and  was 
found  to  be  9.3x10s,  4.2xl06  and  8.4xl07  respectively.  Since 
these  values  are  greater  than  the  external  transition 
criteria  of  5x10s,  TEXSTAN  is  begun  in  the  turbulent  mode. 
The  TEXSTAN  "1  aw-of -the-wal 1"  option  (Crawford  and  Kays, 
1976:25)  is  used  to  speed  calculation.  The  edge  of  this 
region  (y+  =  u+  =  10.)  is  calculated  in  the  INJECT  subroutine 
TLAWDP. 

Results 

The  injection  region  analysis  was  conducted  at  Mach  2.6 
because  of  the  availability  of  the  Fernandez  and  Zukoski  data 
for  5  versus  However,  they  did  not  measure  velocity 
profiles  downstream.  The  following  analysis  is  conducted  at 
Mach  4.38  because  of  the  availability  of  downstream  profiles 
from  Rozycki's  test  program. 

Boundary  layer  development  for  the  uniform  and  '  f' 
function  starting  profiles  are  compared  in  Figure  18.  For 
the  Mach  4.38  conditions,  the  initial  profile  and  the 
profiles  at  3  downstream  locations,  6,  12  and  18  1  heights 
downstream,  are  presented.  The  solid  curves  represent 
profiles  originating  from  the  uniform  profile,  the  dashed 
curves  represent  profiles  originating  from  the  'f'  function. 
For  the  uniform  starting  profile,  the  "law-of-the  wall" 
region  may  be  seen  at  the  bottom  of  the  profile  and  the 
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third-order  polynomial  may  be  seen  at  the  top  joining  the 
uniform  station  e  velocity  to  the  higher  region  d  velocity 
through  the  slip  stream.  In  just  six  boundary  layer  heights 
downstream,  the  uniform  profile  has  already  resolved  itself 
to  the  form  of  the  *  f  *  function  profile  due  to  the  high 
velocity  gradients  at  the  top  and  bottom  of  the  starting 
prof  ile. Thus ,  the  use  of  the  'f'  function  does  not 
significantly  reduce  the  length  of  the  plate  aft  of  the 
injection  region  needed  to  stabilize  the  boundary  layer. 
However,  the  uniform  starting  profile  predicts  greater 
boundary  layer  thicknesses  due  to  the  large  velocity 
gradients  at  the  top  of  the  profile. 

The  boundary  layer  thicknesses  are  more  accurately 
compared  using  the  displacement  thicknesses  which  are 
indicated  on  Figure  18  by  squares.  The  displacement 

thickness  clearly  grows  downstream  and  the  uniform  profile 
displacement  thicknesses  are  greater  than  the  '  f'  profile 
thicknesses  in  each  case. 

In  order  to  validate  the  estimated  3  percent  value  of 
slip  line  thickness,  the  effect  of  slip  line  thickness  on 
boundary  layer  development  was  investigated.  Results  are 
presented  in  Figure  19  for  Mach  4.38  conditions.  The 
thickness  growth  parameter,  o,  is  varied  from  the  baseline  3 
percent  up  to  7  percent.  The  boundary  layer  height 
increases  with  a  since  this  parameter  directly  increases  the 
thickness  of  the  starting  profile;  however,  for  a  130  percent 
increase  in  a,  the  boundary  layer  height  increases  only  6 
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percent.  The  profile  shapes  are  almost  identical.  Thus, 
boundary  layer  development  is  not  sensitive  to  the  assumed 
value  of  slip  line  thickness.  The  original  value  of  3 
percent  is  therefore  used  throughout  this  study. 

Apparently  only  one  set  of  experimental  data  is 
available  for  the  profiles  aft  of  the  injection  region.  In  a 
research  program  (Rozycki),  a  ramjet  inlet  was  forced  to 
ingest  the  missile  forebody  boundary  layer.  Inlet  testing 
was  done  at  the  Air  Force  Academy's  12x12  inch  cross  section 
facility  which  did  not  have  adequate  length  to  reproduce  the 
missile  forebody  boundary  layer.  Boundary  layers  were 
simulated  using  mass  injection  and  profiles  were  measured 
with  pitot  rakes  11  in.  downstream  of  injection  (Figure  20). 
Profiles  for  a  Mach  4.38  test  are  shown  in  Figure  21  and 
Rozycki  has  arbitrarily  laid  power  law  profiles  through  the 
data.  The  circles  represent  the  no  blowing  case.  Two 
blowing  cases  are  plotted,  the  first,  with  5  =  0.627  inches 

99 

is  2.53  times  the  non-blowing  height  and  the  second,  with 
6gg=  0.508  inches,  is  2.02  times  the  non-blowing  height.  The 
data  reasonably  fits  the  power  laws  except  near  the  wall 
where,  contrary  to  intuition,  the  flow  has  higher  velocities 
than  a  non-blowing  boundary  layer  with  the  same  thickness. 
Apparently,  between  station  c  and  e  (this  report's  notation) 
the  flow  rapidly  accelerates  near  the  wall. 

The  velocity  profiles  from  this  study  can  be  compared  to 
Rozycki 's  data  based  on  the  same  99  percent  boundary  layer 
height  but  not  based  on  the  same  injected  mass  flux  because 
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the  author  made  an  error  determining  it.  The  author  had 
available  vendor  information  for  the  porous  plate  mass  flux 
as  a  function  of  Ap  across  the  plate.  For  Ap,  the  author  used 
the  difference  between  plenum  pressure  and  free  stream 
static.  The  actual  pressure  within  the  injection  region  can 
be  typically  double  the  free  stream  value  and  thus  would  have 
yielded  about  half  the  mass  flux. 

Figure  22  compares  the  experimental  data  with  the 
profiles  resulting  from  the  '  f'  and  uniform  starting 
profiles.  The  mass  flow  parameter  X  was  varied  using  INJECT 
until  the  99  percent  boundary  layer  profile  matched  the  data. 
For  6gg=  0.627”,  the  'f'  profile  compares  well  with  the  data 
except  near  the  wall.  For  6gg  =  0.508",  profile  agreement  is 
good  in  the  upper  profile  but  again  falls  off  near  the  wall. 
The  uniform  profile  does  not  represent  the  experimental 
prof i 1 es  wel 1 . 

The  goal  of  this  thesis  is  to  develop  a  general  model 
for  predicting  thicknesses  and  profiles  resulting  from  mass 
addition.  In  the  Introduction,  a  Mach  10,  Nitrogen,  inlet 
test  was  considered  which  required  6  feet  of  extension  plates 
to  generate  a  0.5  in.  boundary  layer  for  use  with  a  1.5  in. 
cowl  height.  The  computational  model  developed  in  this 
report  can  now  be  used  to  model  the  artificial  thickening  of 
the  boundary  layer  to  illustrate  the  use  of  mass  addition. 
The  design  variables  are  the  length  of  the  injection  region, 

1  ,  the  injected  mass  flux  parameter,  X  ,  and  the  length  of 

i  ® 

the  plate  aft  of  injection  to  stabilize  the  profile. 


26 


Nitrogen  is  injected  into  the  Nitrogen  free  stream. 

Retaining  the  overall  plate  length  of  6  ft,  X  is  varied 

00 

until  5  reaches  the  14  degrees  separation  limit.  Lambda  is 
held  at  this  value.  The  injection  length  can  then  be  varied 
to  achieve  various  thicknesses  since  the  starting  boundary 
layer  height  (1  )  is  geometrically  1  -tan6.  Figure  23 

e  i 

presents  boundary  layer  profiles  at  6  ft  for  injection 

lengths  of  0,5,  1.0  and  2.0  ft.  The  0.50  in.  boundary  layer 
height  result  from  the  Introduction  is  plotted  for 
comparison.  Boundary  layer  thicknesses  of  7,  11,  and  17 

times  the  non-blowing  thickness  are  obtained,  respectively. 
One  might  expect  a  1/7  power  law  profile  for  a  non-blowing 

turbulent  boundary  layer  so  these  are  shown  for  the  same  99 

percent  thickness.  The  velocities  in  the  injection  boundary 
layer  are  lower  than  the  1/7  power  law  velocities.  The  shape 
of  the  injection  profiles  does  not  appear  to  change 
downstream.  These  profiles,  having  a  change  in  concavity, 
resemble  boundary  layer  profiles  for  flows  in  a  adverse 

pressure  gradient  (White,  1974:470);  however,  no  pressure 
gradient  exists.  The  corresponding  temperature  profiles  are 
presented  in  Figure  24. 
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V_.  Conclusions 


This  work  has  shown  that  the  phenomena  of  uniform 
injection  into  a  supersonic  stream  can  be  modeled  and 
predicted.  The  phenomena  is  driven  by  a  conservation  of  mass 
in  the  injection  region  and  is  not  sensitive  to  changes  in 
the  free  stream  or  injection  conditions.  The  finite 
thickness  of  the  viscous  contact  surface  and  leading  edge 
curvature  have  the  effect  of  increasing  the  contact  surface 
deflection  angle.  The  injection  region  ('f')  velocity 
profile  predicts  the  boundary  layer  thickness  more  accurately 
than  a  simple  uniform  profile.  The  boundary  layer 
thicknesses  are  not  sensitive  to  the  initial  assumptions  for 
the  slip  line  thickness.  The  predicted  boundary  layer 
profiles  resemble  the  injection  profiles  of  Rozycki ,  however, 
injection  profiles  do  not  resemble  the  1/7  power  law 
turbulent  profile  (A  means  of  rectifying  this  problem  is 
proposed  as  the  second  recommendation). 

The  modeling  of  this  complex  flow  phenomena  has  been 
made  tractable  by  some  assumptions  and  simplifications. 
Sufficient  accuracy  has  been  obtained  however  to  provide  a 
useful  tool  for  the  design  of  boundary  layer  generators. 
This  model  allows  the  designer  to  size  the  system.  Any 
inaccuracies  can  be  compensated  for  during  actual  testing  by 
monitoring  a  boundary  layer  rake  and  adjusting  the  injection 
mass  flow  accordingly. 


VI .  Recommendations 


Two  areas  are  suggested  for  further  study.  The  first  is 
the  effect  of  the  turbulent  mixing  thickness  at  the  contact 
surface  in  the  upstream,  injection  solution.  This  could  be 
analyzed  by  modifying  the  boundary  layer  code  for  two  free 
streams  and  applying  it  to  the  contact  surface  conditions. 
One  problem  would  be  the  starting  conditions.  The  3  percent 
thickness  and  third  order  polynomial  profile  could  be  used  as 
inputs  to  begin  the  solution. 

A  method  of  solving  this  problem  might  be  as  follows. 
First,  INJECT  would  be  run  in  its  current  form  to  provide  the 
conditions  in  the  injection  region  and  region  b.  Region  b 
provides  a  uniform  upper  boundary  condition  for  the  viscous 
code.  The  lower  boundary  condition  would  be  a  function  of  the 
viscous  layer  width  since  y  would  change  as  a  function  of 

up 

x.  Thus,  the  boundary  condition  at  the  lower  edge  of  the 
contact  surface  thickness  would  have  to  be  evaluated  at  each 
space  step.  Once  the  profile  and  thickness  of  the  viscous 
layer  is  determined,  a  displacement  thickness  would  be 
calculated  and  added  to  the  contact  surface  deflection  angle 
in  a  second  run  of  INJECT.  The  process  would  then  be 
repeated  until  little  change  in  the  solution  is  noted. 

The  second  recommendation  for  future  effort  is  to  devise 
a  method  and  computational  model  to  manipulate  the  injection 
boundary  layer  profiles  to  more  accurately  model  the 


non-blowing  1/7  power  law  profiles.  This  might  be 
accomplished  by  inserting  rows  of  holes  aft  of  the  existing 
injection  region  which  injected  high  momentum  flow  angled 
downstream  to  "fatten  up"  the  injection  profile.  Aft  of  this 
row,  an  upstream  facing  row  is  needed  to  slow  the  profile 
near  the  wal 1 . 
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Fig.  2  Inlet  Test  model  with  6  ft  extension  plate 
used  to  create  a  representative  vehicle 
forebody  boundary  layer. 
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Fig.  3  Boundary  Layer  Generator  System  Schematic 


(Hartunian  and  Spancar, 1966: 1306) 
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Fig.  6  Injection  Region  Flowfield  Geometry. 
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Fig.  16  VELOCITY  PROFILE 

Exarr.pie  of  the  *f‘  function 
and  3rd  Order  Polynomial 
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Fig.  17  TEMPERATURE  PROFILE 

Example  of  the  'f'  function 
and  3rd  Order  Polynomial 
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Fig.  18  Velocity  Profiles 

Aft  of  Injection  Region 
Mach=  4.38 
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Fig.  19  Velocity  Profiles 
0.90  n  1  1  inches  Downstream 


co  tn  c- 
o  o  o 


o  o  o 


(ui)a 


49 


500  1000  1500  2000  2500 


Fig.  20  Injection  Hardware  for  Academy  Test, 
Mach=4.38  (Rosycki , 1968 : IV-3 ) 


Pig.  21  Boundary  Profiles  for  Academy  Test, 
Mach=4 . 38  (Rosycki ,1968: IV-16) 
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Fig.  22  Velocity  Profiles  1 1  Inches 
Aft  of  Injection  Region 
M  a=  4.38 
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Appendix  A 

Flowchart  for  INJECT. FOR 
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Appendix  8; 


txxxxxxxxxxxxxxxxxxtxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


x  x 

X  INJECT. FOR  * 

*  * 

*  Capt  Robert  Clausen  * 

*  Air  Force  Institute  of  Technology  * 

*  t 

x  INJECT. FOR  is  a  Fortran  program  which  models  mass  addition  (blowing)  * 

*  on  a  fiat  piate  into  a  supersonic  stream.  This  program  calls  the  boundary  * 

*  layer  code  TEXSTAN  (Kays  and  Crawford)  to  then  calculate  boundary  layer  * 

*  development  aft  of  the  injection  region.  In  order  to  understand  the  * 

*  injection  phenomena  and  analysis  method  of  this  code,  please  obtain  a  * 

*  copy  of  the  AFIT  Thesis  'A  Computational  Model  For  Thickening  Boundary  * 

*  Layers  with  Mass  Addition  for  Hypersonic  Engine  Inlet  Testing".  * 

*  * 


txxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

PROGRAM  INJECT 

Implicit  Real*8  (a-h,m,l,o-z) 

Character*72  filenm, title 
grekpi=  3.141592654 
c 

c —  Open  input  file  'injs.inp'  and  output  file  — 
c 

OPEN  ( 3 .f ile= ’inject . inp ’ ,status=’old’ ) 

Read( 3 ,406  )f ilenm .gamma ,Cp ,R ,  Ma ,pta ,Tta ,Ti ,li .deltid , lambda , 

$  xl.slipfr .title 

406  Format! lx ,/// ,13x ,A14,/,12(  13x,fl0.0,/),/,A72) 

OPEN  (6,file=filenm,status=’unknown’ ) 

Write(6,458)title 
458  Format( A72,//) 

Ta=Tta/(  1 .0+(  gamma-1. 0)/2.0*Ma*Ma) 

pa=  pta/(  1 . 0+(  gamma- 1 . 0  )/2 .  0*Ma*Ma  )**(  gamma/(  gamma- 1.0)) 

deltai=  deltid  *  grekpi/180.0 

mdoti=  lambda*( pa/R/Ta*Ma*DSQRT( gamma*R*Ta )  )*li 


c 

C--  Determine  Mach  wave  angle  (thetam)  as  a  function  of  freestream  Mach  -- 
c  and  use  it  for  the  initial  guess  for  the  iteration, 
c 

thetaMW=  DASIN( 1 .O/Ma ) 
theta«  thetaMU  ♦  10.0*grekpi/l80.0 
c 

c —  Iterate  on  theta  using  a  Newton-Raphson  method  until  mdote*mdoti  — 
c 

iter=0 

dthetas0.02*grekpi/180.0 

10  CALL  MDDIFF(mddifno«,  theta, Ma, pa,  Ti,li,mdoti,deltai, 

$  gamma, Cp,R,  pe.Ue.Te,  pi,Ui,pd,Me,delta,Md) 

IF  (DA8S( mddifnom)  .LE.  (mdoti/1000.0))  GO  TO  11 
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iter=  iter+1 

IF  (iter  .GT.  25)  GO  TO  12 

CALL  MDDIFF(mddifr,  theta+dtheta.Ma.pa,  Ti , li .mdoti .delta! , 

$  gamma, Cp,R,  pe.Ue.Te,  pi ,Ui , pd, Me, delta  ,Md) 

CALL  MD0IFF( mddif 1 ,  theta-dtheta ,Ma  ,pa ,  Ti  ,li  .rndoti  .deltai , 

$  gamma, Cp,R,  pe.Ue.Te,  pi ,Ui ,pd ,Me .delta ,Md ) 
theta=  theta-  (2.0*dtheta*mddifnom)/(mddifr-mddifl ) 

GO  T0  10 

12  Write(6,407)iter-l,mddifnom 
Write!  *  ,407  )iter-i  .mddif  nom 

407  Format!’  The  number  of  iterations  has  exceeded  M2, , 

$  ’  mdiff=’ ,f 10 . 1 ) 

11  Continue 
c 

c—  Write  injection  region  output  — . 

c 

Ua=  Ma*DSQRT(gamma*R*Ta) 

Write!  6, 503  )Ma,Ua,pa,pta,Ta,Tta 

503  Format!  ’  Ma=  ’  ,f 5 .2,2x ,  'Ua( ft/s )=’  ,f6.1 ,2x ’pa( psf  )=’  ,f6.1 ,2x , 

$  ’pta( psf  )=’  ,fl0.1,2x,/,’  Ta(R  )=  ’ ,f6.1 ,2x , ’Tta( R)=  ’  ,f6.1 ,/) 

Write!  6,502  )delta*57.296 

502  Format!’  Contact  surface  ('wedge')  angle  delta! deg )=’ ,f5 .1 ,/ ) 

Write! 6,504  )mdoti .lambda ,pi ,Ui  ,Ti .deltid.li 

504  Format!  ’  rndoti! slug/s )=’  ,f7. 5, 2x, ’lambda*’  ,f6. 4, 2x,  ’pi( psf  )=  ’ , 

$  f 6 . 1 ,2x ,  ’Ui( ft/s )=  ’  ,f  6 . 1 ,2x ,/ ,  ’Ti!  R  )=  \f  6 . 1 ,2x , 

$  ’  deltai!  deg  )= ’f5.  l,2x,’li(  ft  )=’  ,f6.4,/) 

Td=  T  t  a  /( 1 . 0+( gamma- 1 . 0 )/2 . 0*Md*Md ) 

Ud=  Md*0SQRT( gamma*R*Td ) 

ptd=  pd*(  Tta/Td  )**(  gamma/!  gamma- 1.0)) 

wr i te( 6 , 505  )Md ,Ud,pd,ptd,Td,Tta 

505  Format!  ’  Md=’  ,f 5 .2 ,2x ,  ’Ud(ft/s)=’  ,f6.1 ,2x,  ’pd(psf  )=’  ,f6.1 ,2x , 

$  ’ptd( psf  )*’,f8.1,2x,/,’  Td(R)=’,f6.1,2x,’Ttd(R)=’,f6.1,/) 

le=  li*0TAN(  delta) 

Write!  6 ,506  )Me  ,Ue  ,pe  ,Te ,  le 

506  Format!  ’  Me=’,f5.2,2x,’Ue(ft/s)=’,f6.1,2x,’pe(psf  )»’,f6.1,2x, 

$  ’Te(R)*’  ,f6.1 ,2x ,/,  ’le(  ft  )=’  ,f6.4,/) 

c 

c —  Calculate  the  Reynolds  number  at  station  e  - 

c 

Re=  ud*(li+le)/visc(Td,pd) 

Write!  6, 507  )Re 

507  Format!  ’  Reynolds  Nu.  at  Station  E  =  ’,el0.4,/) 
c 

c —  Calculate  the  boundary  layer  thickness  at  these  conditions  on  a  plate 
c  of  the  same  length  for  comparison, 
c 

del99=  0.37  *  xl  *  (1.0d0/(Ua*xl/visc(Ta,pa)))**.2 

deli»  del99/8.0 

del2=  del 99  *7.0/72.0 

Write!  6,508  )xl  ,del99*l2 .  ,dell*12 .  ,del2*12 . 

508  Format!’  8oundary  layer  thicknesses  at  xl(ft)*\f5.2,’  without 
linjection’,/, ’  delta99( in)»’,f7.4,4x, 

$  ’displacement  deltai!  in  )«’,f7. 4, 4x, ’momentum  delta2(in)=’,f7.4,/) 
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c 

c—  Determine  fup  from  slipstream  thickness  — 
c 

c=  grekpi/2. O/delta/delta 

IF(0SIN( delta)  .GT.  slipfr/2.0)  GO  TO  26 

Ur ite( * ,421 )  delta*57.3 

421  Format!’  delta=  ’ ,f5.2, ’which  causes  the  slipstream  to  occupy’, 

$  ’  the  entire  region  ”e’.’,/, 

$  ’  Increase  the  injected  mass  flow  (lambda).’) 

STOP 

26  fup*  slipfr/2.0/DSIN( delta) 
hti=  Cp*Ti  +  Ui**2/2 .0 

c 

c—  Find  the  *f ’  profile  which  yields  conservation  of  mass  at  station  e.  --- 
c  Iterate  on  u  at  fup  (uup)  until  the  massflow  calculated  from 

c  integrating  the  'f'  profile  and  slipstream  profile  equals  mdote. 

c  A  Newton-Raphson  iterative  skeme  is  used, 
c 

duup=  0.0001*ud 
uup=  0.84*ud 

27  Teup=  (hti  -  uup«2/2.0)/Cp 
IF  (Teup  .GE.  0.0)  Go  To  28 
uup=  0.89*uup 

Go  To  27 

28  CALL  MDEDFF(mdenom,uup,fup,ud,Td,c,le,hti,mdoti,Cp,R,pe, 

*  au,bu,cu,du,  aT,bT,cT,dT,  Vw) 

IF( DABS( mdenom)  .LE.  (mdoti/1000.0))G0  TO  31 

CALL  MDEDFF( mder ,uup+duup ,fup ,ud  ,Td ,c , le  ,hti .rndoti ,Cp ,R ,pe , 

$  au.bu.cu.du,  aT,bT,cT,dT,  Vw) 

CALL  MDEDFF( mdei ,uup-duup .fup ,ud  ,Td  ,c , le ,hti .mdoti ,Cp ,R ,pe , 

$  au,bu,cu,du,  aT,bT,cT,dT,  Vw) 
uup=  uup  -  (2.0*duup*mdenom)/(mder-mdel) 

GO  TO  27 
c 

c—  Oetermine  y  at  the  edge  of  the  law-of-the-wall  region  so  that  this  option 
c  can  be  used  in  TEXSTAN.  A  Newton-Raphson  iterative  method  is  used. 

c 

31  yup=  le*(1.0-fup) 
yd=  le*(1.0+fup) 
ylawwa*  0.05*yup 
dylaw*  0.0001*yup 
viscw=visc(Ti,pe) 

43  ylawnom*  YLAW0f(  ylawwa,  Vw,le,c,viscw) 

IF( DABS(ylawno«)  .LE.  0.1)  GO  TO  44 
ylawwa*  ylawwe-(2.0*dylaw*ylawnom)/ 

$  ( YLAW0F( ylawwe+dylaw.Vw.le ,c ,viscw )- 
$  YLAUOF( y lawwa-dy law ,Vw , le ,c ,viscw ) ) 

GO  TO  43 
c 

c—  Write  the  velocity  profiles  to  file  ’PROFILE’  for  input  into  Texstan  — 
c  y  in  ft,  Velocity  in  ft/s,  Static  entha  *  in  BTU/lbm. 
c 

44  OPEK  5,file= ’PROFILE ’ .status* ’unknown’ ) 


58 


*rite(5,479  )titie 
479  rormat( A64 ) 
y0=  O.OdO 
.0=  O.CdO 

Ur ite(  5,427  )yO ,uO ,Cp*Ti/778 .0/32 .2 
427  Format(3(el0.4)) 
c 

c--  ’f  ’  profile  region  -- 

c 

Nlower=  40 

dely=(  yup-y lawwa  )/D8LE( Niower ) 

Oo  152  i=l, Niower 
a  i= i 

ye=ai*dely+ylawwa 
fe=i  ,0-ye/le 

uep=  Vw*D$QRT( 1 ,0-fe**2-2.0*c*DL0G(  fe)) 

ne=  ( hti-Vw**2*(  1 .0-fe**2-2.0*c*0LOG(  fe )  )/2 .0  )/778  ./32 .2 

152  klrite(5,427)ye,uep,he 

r 

c--  slipstream  region  — 

c 

Nsi ip=  30 

dely=( yd-yup  )/DBLE(  Nsl ip ) 

Oo  153  i=l,Nslip 
ai=  i 

ye=  ai*dely+yup 

uep=  au*ye**3  +  bu*ye**2  +  cu*ye  +  du 

he=  (aT*ye**3  +  bT*ye«*2  +  cT*ye  +  dT)*Cp/778./32.2 

153  Ur ite( 5 ,427 )  ye.uep.he 
N=  Niower  +  Nslip 
Rewind  5 

Close!  5) 
c 

c—  Determine  the  momentum  and  energy  of  the  fluid  at  station  ’e’  from 
c  the  1-0  properties  and  compare  them  to  the  momentum  and  enery 
c  calculated  by  numerically  integrating  the  'V  velocity  profile, 

c  This  serves  as  a  check  to  determine  how  closely  the  assumed  *f  * 

c  profile  conserves  mass,  momentum  and  energy.  (Note  that  the 
c  assumed  ’f '  profile  at  V  was  defined  by  conserving  massflow). 
c 

c—  Calculate  the  momentum  and  enery  from  1-0  properties 

c 

momlD=  pe/R/Te*Ue**2*le  +  pe*le 
E1D3  pe/R/Te*Ue*le*(Cp*Te  +  Ue**2/2.0) 
c 

c—  integrate  the  momentum  and  energy  profile  in  region  e 
c—  integrate  the  region  below  the  slipline 

c 

df=  (1.0-fup)/200.0 
suml3  0.0 
sum8=  0.0 
Oo  101  k=l  ,200 
ak=  k 
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*=  ( aK-C  .5  )*df  +  fup 

ul=  Vw*D$QRT( 1 .0  -  f**2  -2.0*c*0L0G(f )) 

suml=  suml  +  ul**2/(hti-ul**2/2.0)*df 

101  sum8=  sum8  +  ul  *(1.0  +  ul**2/2.0/(hti-ul**2/2.0))*df 
moml=  pe*le*Cp/R*suml 

Ei=  pe*le*Cp/R*sum8 

c—  irtegrate  the  massflow  in  the  slipline - 

c 

dy=  ( le-yup)/100.0 
sum2=  0.0 
sum9=  0.0 
Do  102  k=l ,100 
ak=  k 

y=  (ak-0.5)*dy  +  yup 
us=  au*y**3  +  bu*y**2  t  cu*y  +  du 
Tes=  aT*y**3  +  bT*y**2  +  cT*y  +  dl 
sum2=  sum2  +  us**2/Tes*dy 

102  sum9=  sum9  +  us*(Cp  +  us**2/2.0/Tes)*dy 
moms1  pe/R*sum2 

Es=  pe/R*sum9 

Horae*  moml  +  moms  +  pe*le 
Ee=  El  +  Es 

perdfm=  DA8S( momlD-mome )/momlD 
perdfE=  DA8S(E10-Ee)/E10 

klr ite( 6 ,466)  momlD ,mome ,perdf m*100 . ,£10/778 . ,Ee/778 . ,perdfE*100 . 

466  Format(  ’  Momentum  at  Station  V’,/,’  1-0  (lbf  )*’,f8.2, 

$  5x , ’  "f*  Profile  (lbf)=’,f8.2,5x,/,’  percent  difference  =’,f7.1, 

$  //,’  Energy  at  Station  *e'\/,’  1-0  (Btu/s)=\f8.2,5x, 

$  ”f“  Profile  (Btu/s)=’ ,f8.2,5x,/,  ’  percent  difference^ ,f7.i ,//) 
c 

c—  Set  up  the  TEXSTAN  input  file - - 

c 

xu=  li+le 

hw=  Ti*Cp/778./32.2 
hd=  Td*Cp/778./32.2 

CALL  INPUTD(  title, xu,xl,pd,ud,hw,hd,N) 

c 

c—  Start  TEXSTAN - 

c 

CALL  TEXSTAN 

STOP 

END 

c 

Cttttttt ttt**t*t****ttttt*t*tttt*tt*tttt*t*tt*tt*t*ttttt*tt*t**ttttttt*t** 

c 

SUBROUTINE  M00IFF(mddif ,  theta, Ma, pa,  Ti ,li .mdoti .deltai , 

*  gamma, Cp,R,  pe.Ue.Te,  pi ,Ui ,pd ,Me .delta ,Hd ) 

Implicit  Real*8  (a-h,m,l,o-z) 
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ceita=  OATAN(  1 .0/(  ( gamma+1 . 0  )*Ma*Ma/2 . 0/( ( Ma*DSIN( theta )  )**2-l .0 ) 

$  -1.0  )/OTAN( theta ) ) 
lc=  li*DTAN(  delta) 

?b=  pa*(2.O*gamma*Ma*Ma*0SIN(theta  )**2  -  ( gamma-1. 0))/(  gamma+1. 0) 
pi=  pb 

Ui=  mdoti*R*Ti/pi/li 
c 

c —  Solve  *'or  the  conditions  in  region  ‘d‘  after  the  expansion. - 

c 

MaSIN2=  ( Ma*DSIN(  tneta )  )«2 

Mb=  DSQRT((( gamma+1 ,0)**2*Ma*Ma*MaSIN2  -  4.0*(Masin2-1.0) 

$  *( gamma*Masin2+l .0 ) )/( 2 .0*gamma*Masin2-(  gamma-1 .0 ) )/ 

$  ((gamma-1 ,0)*Masin2+2.0)) 

ptp=  ?b*(  1 .C+(  gamma-1 .0  )/2 .0*Mb*Mb )**( gamma/(  gamma-1 .0 ) ) 

CALL  EXPANS( Mb, ptb, delta,  Md,pd,  gamma) 

c 

c —  Solve  the  lower  control  volume  for  pe,  Ue,  and  Te.— 
c 

le=  lc 
pe=  pd 

Ue=  (2.0*mdoti*(Cp*Ti  +  Ui*Ui/2.0))/ 

$  (pb*lc  -  pe*le  +  mdoti*Ui*C0S(deitai )  +  2.0*pe*le*Cp/R) 

Te=  (pe*Ue**3*le)/2.0/R/(mdoti*(Cp*Ti+Ui*Ui/2.0)  -  pe*Ue*le*Cp/R) 

Me=  Ue/DSQRT(gamma*R*Te) 
mdote=  pe/R/Te*Ue*le 
mddif=  mdoti  -  mdote 
RETURN 
END 
c 

QtXXXXXXXXXXXXXXXXXXXXXXXXXXXXXtXXXtXttXtttttttXXXXtXXXXXXXXXXXXXXXXXXXtXXXX 

c 

SUBROUTINE  EXPANS  (Ml.pt.d,  M2,p2,  gamma) 

Implicit  Real*8  (a-h,m,l,o-z) 

c 

c  This  subroutine  calculates  the  Mach  number  and  pressure  after  a 
c  Prandtl-Meyer  expansion  of  angle  *d‘.  Function  subroutines  "DIF’ 
c  and  "NU’  are  part  of  this  subroutine, 
c 

dM=  0.05 
M2=  Ml  +.2 

10  DIFN0M=0IF(  gamma, Ml  ,M2,d) 

IF(DABS(DIFNOM).LT. 0.0017)  GO  TO  100 
M2=  M2-( DIFNOM*2.0*dM )/( DIF{  gamma ,M1 ,M2+dM ,d )- 
$  DIF( gamma  ,M1 ,M2-dM  ,d ) ) 

GO  TO  10 

100  p2=  pt/(  1 .0  +  ( gamma-1 ,0)/2.0*M2**2)**(gamma/( gamma-1.0)) 

RETURN 

END 

C 

FUNCTION  DIF( gamma, Ml ,M2,d) 

Implicit  Real*8  (a-h,m,l,o-z) 

Rea 1*8  NU 
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D  IF=NU(  gamma  ,M2 )-( d+NU(  gamma  ,M1 ) ) 

RETURN 

END 

C 

FUNCTION  NU(  gamma  ,M) 

Implicit  Real*8  (a-h,m,l,o-z) 

Rea  1*8  NU 

b=(  gamma+l)/(  gamma-1) 
c=M**2-1.0 

NU=(  D$GRT(  b  )*DATAN(  DSQRT(  c/b )  )-DATAN(  DSQRT(  C ) ) ) 

RETURN 

END 

c 

c  *********************************************************************** 

c 

SUBROUTINE  MDEDFF( mdedif ,uup ,fup ,ud ,Td ,c , le ,hti .rndoti ,Cp ,R ,pe , 

$  au,bu,cu,du,  aT,bT,cT,dT,  Vw) 

Implicit  Real*8  (a-h.m.l.o-z) 

c 

c  This  subroutine  numerically  integrates  the  mass  flow  in  the  ’f’  profile 
c  and  slipline  profile  region, 
c 

yup=le*(  1.0-fup) 
yd=  ie*(1.0+fup) 

Vw=  DSQRT(  uup**2/(  1 .0-fup«2-2 .0*C*0L0G(  f  up ) ) ) 

Teup=  (hti-uup**2/2.0)/Cp 
zT=  -Vw**2/Cp/le*(fup+c/fup) 

CALL  COEF(yd,yup,zT,Td,Teup,  aT,bT,cT,dT) 

zu=  Vw*( fup+c/fup  )/le/DSQRT( 1 .0-fup**2-2 .0*c*DL0G( fup ) ) 

CALL  COEF(yd,yup,zu,ud,uup,  au,bu,cu,du) 
c 

c —  integrate  the  mass  flow  profile  in  region  e  - 

c 

c—  integrate  the  region  below  the  slipiine - 

c 

df=  (1.0-fup)/200.0 
suml=  0.0 
Do  101  k=l  ,200 
ak=  k 

f=  (ak-0.5)*df  +  fup 

sub=  1.0  -  f**2  -2.0*c*0L0G(f ) 

101  su«l=  suiil  +  Vw*DSORT(sub)/(hti-Vw**2/2.0*sub)*df 
mdotl*  pe*le*Cp/R*suml 

c 

c—  integrate  the  massflow  in  the  si  ip  line - 

c 

dy=  ( le-yup)/100.0 
sum2-  0.0 
Do  102  k=l  ,100 
ak=  k 

y=  (ak-0.5)*dy  +  yup 

us*  au*y**3  t  bu*y**2  +  cu*y  +  du 

Tes*  aT*y**3  +  bT*y**2  +  cT*y  +  dT 
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102  3um2=  sum2  +  us/Tes^dy 
■naots=  ?e/R*sum2 
,ndote=  indot  1  +  mdots 
mdedif =  mdoti-mdote 
RETURN 
END 
c 

zxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxtxxxxxxxxxxxxxxxxxxxxtxxxxxxxxxxxxxxx 

r 

SUBROUTINE  CCEF(yd,yu,b2,b3,b4,  aa,bb,cc,dd) 

Implicit  Reai*8  (a-h,m,l  ,o-z) 

Real*8  A(4,4),b(4) 
c 

c  ’"'.is  subroutine  inverts  a  four  by  rour  matrix  in  order  to  solve  for 
c  the  coeficients  of  the  third  order  polynomial  used  to  model  the 

c  velocity  and  temperature  profile  through  the  slipline. 

c 

A(l,l)=  3.*yd**2 
A(  1 , 2  )=  2.*yd 
A(1.3)=  1.0 
A(  1 ,4  )=  0.0 
A(  2 ,1  )=  3.*yu**2 
A(  2 ,2  )=  2.*yu 
A(  2 ,3  )=  1.0 
A(  2 ,4  )=  0.0 
A(  3 ,1  )=  yd«3 
A(3,2)=  yd«2 
A(  3,3  )=  yd 
A(3,4)=  1.0 
A(  4 , 1  )=  yu**3 
A(  4 ,2  )=  yu«2 
A(  4 ,3  )=  yu 
A(  4 ,4  )=  1.0 
b(  1  )=  0.0 
b(2)=  b2 
b(3)=  b3 
b(4)=  b4 

Do  122  k=l ,3 
Do  122  i=k*l,4 
amult=  A(i,k)/A(k,k) 

Do  123  j=k+l  ,4 

123  A(i,j)=  A(i,j)-A(k,j)*amult 
A(i,k)=0.0 

122  b(i)=  b(i)-b(k)*amult 
c 

c -  back  solve  - 

c 

b(  4  )=b(  4  )/A(  4,4) 

Do  124  k=3,l ,-l 
sum=0.0d0 
Do  125  j=k+l ,4 

125  sum=  A(k,j)*b(  j)+sum 
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124  b(x)=  ( b(  k  )-sum  )/A(  k ,  k ) 

dd=b(  4 ) 
cc=b(3) 
ob=b(  2 ) 
aa=b(l) 

RETURN 

END 

C 

ctxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c 

FUNCTION  YLAWDF( y ,Vw,Ie ,c .vise ) 

Implicit  Real*8  (a-h.m.l.o-z) 

c 

c  This  subroutine  returns  the  difference  between  the  RHS  and  LHS  of 
c  the  law-of-the-wali  equation  u+*y+  =  U*y/nu  to  the  portion  of 

c  the  main  program  seeking  to  find  the  value  of  y  which  makes  this 

c  difference  zero. 

c 

ypius=  10. 
uplus=  10. 

IF( ( y  .IT.  O.OdO)  .OR.  (y  .GT.  ie))G0  TO  83 
f=  1.0-y/le 

u=  Vw*DSQRT(  1 .0-f  «2-2 .0*c*DL0G(  f ) ) 
ylawdf=  u*y/visc  -  uplus*yplus 
Return 

83  Urite(*,409)y , le 

409  Format( '  y  is  out  of  bounds  in  Subroutine  YLAUDF.  y  must  be  0.0  ( 

$  y  <  le  y=  ’ ,dl2.5, ’  le»  ’,dl2.5) 

STOP 

End 

c 

cxzttttttttttttttttttttttxtttttttttttttxxxttttttttttttttxttttttttttttttttt 

c 

SUBROUTINE  INPUTD( title , xu , xl , pd ,ud ,hw , hd , n ) 

Implicit  Real*8(a-h,o-z) 

INTEGER  GEOM .FLUID .SPACE .OUTPUT .MODE 
Character*72  title 
c 

c  This  subroutine  sets  up  the  TEXSTAN  input  file  'DATA'  which  governs 
c  geometry,  boundary  conditions,  flow  conditions,  turbulence  models, 

c  and  input/output  instructions.  See  TEXSTAN  manual  for  details, 

c 

0PEN( 7 ,f ile= ’DATA’ ,status= ’unknown ’ ) 

WRITE  (7,100)  title 
100  FormatC  ’ ,A72,’  ’) 

geom=l 

mode=2 

fluid=2 

neq=2 

kex=2 

kin=2 
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kin=l 

kent=l 

kstart=0 

<kd=10 

WRITE  (7,110)  GEOM, MODE, FLUID  ,NEQ,N,KEX, KIN, KENT .KSTART.KKD 

deltax=0.4 

retran=200.0 

£ra=0.05 

enfra=0.005 

gv=0.0d0 

WRITE  (7,120)  XU  ,XL  .OELTAX  .RETRAN  ,FRA  .ENFRft  ,GV 

X  *  "  * 

write(7,no)il,il 


rhoc=l  .0 

viscoc=.10d-04 

prc=0.7 

JRITE  (7,120)  Pd,RH0C,VISC0C,PRC 


12=2 

Write(7,ll0)i2,il 

Write(7,110)i2,il 

al=1.0d0 

aO=O.OdO 

Ur  ite(  7 , 120  )xu  ,al  ,aO  ,aO  ,aO 
Urite(7,120)xl,al,a0,a0,a0 
Write(7,120  )aO,aO,hw 
Urite(7,120)ud,a0,hd 
Urite(7,120)a0,a0,hw 
Ur  ite(  7 , 120  jud.aO.hd 

ak=  0.41 

almgg=0.085 

fr=0.01 

aq=0  .OdO 

bq=O.OdO 

ypMX*50.0 

yp«in*20.0 

WRITE  (7,120)  AK , ALMGG ,  FR ,  AQ ,  80 , YPMAX , YPMIN 

apl 1=25.0 

bpll=0.0 

signal=1.0 

WRITE  (7,120)  APLL.BPLL, SIGNAL 

pplag=4000.0 

prt=0.90 

Write(7,120)ppLag,prt 

cg=32.2 
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c  >773.0 

wr  ite(  7 , 120  )cg  ,c  j  ,a0  ,a0  ,a0  ,a0  ,a0 
numrun=l 
space= 10000 
output=2 
kl=0 
.<2=3 
<3=0 
k4=0 
<5=0 
<6=0 
k7=0 
k8=0 
k9=99 
<10=0 
kll=0 
k!2=0 
k!3=0 

WRITE( 7,110  )NUMRUN .SPACE  .OUTPUT  ,K1  ,K2  ,K3  ,K4  ,K5 ,K6 ,K7 ,K8 ,K9  ,K10 , 

$  K11.K12.K13 

110  Format!  16(  15)) 

120  Format! 8( G10 .4 ) ) 

Rewind  7 
Close  (7) 

Return 

End 

c 

C 

FUNCTION  VISC(TR.p) 

Implicit  Real*8  (a-h.l.ra.o-z) 

c  Calculation  of  kinematic  viscosity.  Dynamic  viscosity  (mu)  is  determined 
c  as  a  function  of  temperature  below  3000  K  from  the  relation  in 

c  Bertin  and  Smith.  Nu  is  not  dependent  on  pressure  in  this  range. 

c 

R=  1717.93 
TK=  TR/1.8 

mu=  1.458d-06  /  47.88026  *  tk«l  .5/(  tk+110.4 ) 

visc»  mu/p*R*TR 

Return 

End 
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Appendix  C 


Sample  Input 


Input  data  for  “inj.for" 


output  £ile=  Rozycki. out 

1 . 

gamma =  1.4 

2. 

Cp=  6008.61 

f t*ibf /slug/R 

3. 

R=  1717.9344 

f t*lbf /slug/R 

4. 

Ma=  4.38 

5. 

Pta=  29088. 

psf 

6. 

"ta=  560.0 

R 

7. 

Ti=  560.0 

R 

9. 

li«  0.25 

ft 

9. 

deltid=  90.0 

deg 

10. 

lambda=  0.0153 

11. 

xl=  0.9405 

ft 

12. 

sl ipf r=  0.03 

13.  TEXSTAN  output  file  title  (72  char )= 
Rozycki : Academy  Test 
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Appendix  9 


Saaple  Output 


Rozycki : Acadeay  Test 


la-  4.38  Ua(ft/s )=2311.3  pa(psf)?  116.9  pta(psf)*  29088.0 
Ta(R)«  115.8  Tta( R )=  560.0 

Contact  surface  ('wedge')  angle  delta( deg )=  8,1 

adoti( slug/s )=0 . 00520  laabda=0 .0153  pi(psf)*  265.6  Ui( ft/s )*  75.3 
Ti(R)®  560.0  deltai( deg )=  90.0  li( f t )=0 . 2500 

Md=  4.34  Ud( ft/s )=2306 .4  pd(psf)-  117.0  ptd(psf)*  27516.1 
Td(R)»  117.7  Ttd(R)*  560.0 

He=  0.95  Ue( ft/s ;=1017 .5  pe(psf)*  117.0  Te(R)*  474.3 
le( ft  )=0 .0356 

Reynolds  Nu.  at  Station  E  =0 .4163E+07 

Boundary  layer  thicknesses  at  xl( ft )=  0.94  without  injection 

delta99( in )=  0.1550  displacement  deltal( in )-  0.0194  aoaentue  delta2( in)=  0.0151 

Momentum  at  Station  V 

1-0  ( lbf )=  9.45  T  Profile  ( lbf )■  10.52 

percent  difference  =  11.4 

Energy  at  Station  V 

1-0  (Btu/s)*  22.49  *f*  Profile  (Btu/s)*  22.44 

percent  difference*  0.2 


1  Rozycki: Academy  Test 

GEOMETRY  MODE  FLUID  NEQ  N  KEX  KIN  KENT  KSTART  KKD 
1  2  6  2  70  2  1  1  0  10 

ENTRAINMENT  8ASE0  ON  BEHAVIOR  OF  ALL  EQUATIONS. 

XU  XL  OELTAX  TRAN.  RE  NO.  FRA  ENFRA  GRAVITY  CONST. 

0 . 2856E+00  0 .9405E+00  0.400QE*00  0.2000E+03  0.5000E-01  0.5000E-02  O.OOOOE+OO 

300Y-F0RCE  SOURCE) 1 )  SOURCE) 2)  SOURCE) 3)  SOURCE) 4)  SOURCE) 5) 

1  1 
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PRESSURE  DENSITY  VISCOSITY  PftC(l)  PRC{ 2 )  PRC( 3 ) 

'HE  “L JIO  IS  NITROGEN 

0.11 70E+03  0.1000E+01  0.1000E-04  C  .7000E+00 

30UNDARY  CONDITIONS  ALONG  I-  ANO  E-SURFACES 
NXSC  "YP8C1  'YP8C2  TYP8C3  TYP8C4  TYP8C5 
2  1 

2  1 

H  X(H)  RU(H)  AUXl(H)  AUX2(M)  AUX3(N)  06(H)  AH(K.J)  FJ(K.l.H)  FJ(K,2,H 

1  0 . 28S6E+00  0.1000E+01  0 .OOOOE+OO  0 .OOOOE+OO  0. OOOOE+OO  0 .2306E+04  O.OOQOE+OQ  0 .1343E+03 

0 . OOOOE+OO  0.2822E+02 

2  0 .94C5E+00  0.1000E+01  O.OOOOE+OO  O.OOOOE+OO  O.OOOOE+OO  0.2306E+04  O.OOOOE+OO  0.1343E+03 

O.OOOOE+OO  0.2822E+02 


INITIAL  PROFILES 

I  Y(I)  0(1)  F(  1 ,1 )  F(  2 , 1 )  F(3,I) 

I  O.OOOOOOE+OO  O.OOOOOOE+OO  0.134300E+03 

3  0.228200E-02  0.262800E+03  0.133100E+03 

4  0.303900E-02  0.304900E+03  0.132600E+03 

5  0.379500E-02  0.342600E+03  0.132100E+03 

6  0.455100E-02  0.377400E+03  0.131600E+03 

7  0.530700E-02  0.410000E+03  0.131100E+03 

8  0.606300E-02  0.440800E+03  0.130500E+03 

9  0.681900E-02  0.470400E+03  0.130000E+03 
10  0.757500E-02  0.498900E+03  0.129500E+03 

II  0.833200E-02  0.526600E+03  0.128900E+03 

12  0.908800E-02  0.553600E+03  0.128300E+03 

13  0.984400E-02  0.580100E+03  0.127700E+03 

14  0.106000E-01  0.606200E+03  0.127100E+03 

15  0.113600E-01  0.631900E+03  0.126500E+03 

16  0 . 121 100E-01  0.657300E+03  0.125800E+03 

17  0.128700E-01  0.682600E+03  0.125100E+03 

18  0.136200E-01  0.707800E+03  0.124400E+03 

19  0.143800E-01  0.733000E+03  0.123700E+03 

20  0.151400E-01  0.75820QE+03  0.123000E+03 

21  0.158900E-01  0.783500E+03  0.122200E+03 

22  0.166500E-01  0.808900E+03  0.121400E+03 

23  0.174100E-01  0.834500E+03  0.120500E+03 

24  0.181600E-01  0.860400E+03  0.119700E+03 

25  0.189200E-01  0.886700E+03  0.118700E+03 

26  0.196700E-01  0.913300E+03  0.117800E+03 

27  0.204300E-01  0.940400E+03  0.116800E+03 

28  0 . 21 1900E-01  0.968100E+03  0.115700E+03 

29  0.219400E-01  0.996500E+03  0.U4600E+03 

30  0.227000E-01  0.102600E+04  0.113400E+03 

31  0.234500E-01  0.105600E+04  0.112200E+03 

32  0.242100E-01  0.108700E+04  0.U0900E+03 

33  0.249700E-01  0.111900E+04  0.109400E+03 

34  0.257200E-01  0.115200E+04  0.107900E+03 

35  0.264800E-01  0.U8800E+04  0.106300E+03 


)  FJ(K,3,N) 
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36  0.272400E-01  0 .122500E+04  0.104500E+03 

37  0.27990QE-01  0.126400E+04  0.102500E+03 

38  0.287500E-01  0.130600E+04  0.100400E+03 

39  0.295000E-01  0.135200E+04  0.979600E+02 

40  0.302600E-01  0.140200E+04  0.952200E+02 

41  0.310200E-01  0 . 1 45700E+04  0.920600E+02 

42  0.317700E-01  0 .152000E+04  0.883200E+02 

43  0.320200E-01  0.154400E+04  0.868400E+02 

44  0.322800E-01  C.156900E+04  0.851700E+02 

45  0.325300E-01  0.159700E+04  0.833100E+02 

46  0.327800E-01  0.162600E+04  0 .812900E+02 

47  0 . 330300E-01  0.165700E+04  0.791100E+02 

48  0.332900E-01  0 .168800E+04  0.768100E+02 

49  0.335400E-01  0.172100E+04  0.744000E+02 

50  0.337900E-01  0.175500E+04  0.718900E+02 

51  0.340400E-01  0.178900E+04  0.693000E+02 

52  0.343000E-01  0 . 182400E+04  0.666500E+02 

53  0.345500E-01  0.185900E+04  0.639600E+02 

54  0.348000E-01  0.189400E+04  0.612500E+02 

55  0.350600E-01  0.192900E+04  0.585300E+02 

56  0.353100E-01  0.196400E+04  0.558200E+02 

57  0.355600E-01  0.199800E+04  0.531400E+02 

58  0.358100E-01  0.203100E+04  0.505000E+02 

59  0.360700E-01  0.206300E+04  0.479300E+02 

60  0.363200E-01  0.209500E+04  0.454300E+02 

61  0.365700E-01  0.212400E+04  0.430400E+02 

62  0.368200E-01  0.215300E+04  0.407600E+02 

63  0.370800E-01  0.217900E+04  0.386200E+02 

64  0.373300E-01  0.220400E+04  0.366200E+02 

65  0.375800E-01  0.222600E+04  0.348000E+02 

66  0.378300E-01  0.224600E+04  0.331600E+02 

67  0.380900E-01  0.226400E+04  0.317200E+02 

68  0.383400E-01  0.227900E+04  0.305100E+Q2 

69  0.385900E-01  0.229000E+04  0.295300E+02 

70  0.388400E-01  0.229900E+04  0.288200E+02 

71  0.391000E-01  0.230500E+04  0.283700E+02 
73  0.393500E-01  0.230600E+04  0.282200Et02 

TURBULENCE  CONSTANTS 


KAPPA 

LAH80A 

FR 

AQ 

BO 

MAX  YPL  WF 

MIN  YPL  WF 

0.4100E+00 

0 .8500E-01 

0.1000E-01 

0 . OOOOE+OO 

0.  OOOOE+OO 

0.5000E+02 

0.2000E+02 

APL 

XXX 

SIGNAL 

PPLAG 

PRT(l) 

PRT(2) 

PRT( 3 ) 

0.2500E+02 

0. OOOOE+OO 

O.iOOOE+Oi 

0.4000E+04 

0.9000E+00 

DIMENSIONING 

SYSTEM  CONSTANTS 

ARBITRARY  CONSTANTS 

G-SU8-C 

J-SUB-C 

AXX 

exx 

CXX 

DXX 

EXX 

0.3220E+02 

0 . 7780E+03 

0. OOOOE+OO 

0. OOOOE+OO 

0.  OOOOE+OO 

0. OOOOE+OO 

0. OOOOE+OO 

NO.  OF  RUNS 

OF  DATA 

PRINTOUT  SPACING  OUTPUT  OPTION 

K1  K2 

K3  K4  K5  K6  K7 

K8  K9 

K 10 

[  K12  K13 

1 

tut 

2 

0  3 

0  0  0  0  0 

0  99 

0 

0  0  0 
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IF  REM  IS  ,£SS  THAN  ABOUT  6000,  LANBOA  IS  MODIFIED  BY  AN  INTERNAL  CORRELATION 
PRT  (FOR  THE  ENERGY  EQ)  NEAR  A  WALL  IS  3EING  EVALUATED  BY  AN  INTERNAL  CORRELATION 

ROUTINE  LANSU8  HAS  BEEN  CALLEO 
N  HAS  SHIFTED  TO  69  AT  INTG  -  0 


ROUTINE  .ANSU8  HAS  3EEN  CALLED 
N  HAS  SHIFTED  T0  68  AT  INTG  =  0 


ROUTINE  lANSU8  HAS  BEEN  CALLED 
N  HAS  SHIFTED  TO  67  AT  INTG  =  0 

FLOW  IS  TURBULENT 

USING  PRANDTL  MIXING-LENGTH  TURB .  MODEL 

USING  VAN-DRIEST  EXPONENTIAL  DAMPING  FUNCTION  IN  THE  INNER  REGION 
MODIFYING  LAMBDA  AT  LOU-RE 


INTG  XU  UGU  X  F  REH  CF2  H  REH  ST  F(l,MALL)  A 

ME  EXX 

*  0  0 . 2856E+00  0.2306E+04  -.367E-09  O.OOOE+OO  0.5494E+04  0.7169E-04  0.3207E+01  0.6040E+05  0.1461E-02  0.1343E+03  0.00 

OE+OO  0 . 223E+03 


I 

Y(I) 

U(I) 

F(l.I) 

F(2,I) 

F(  3 , 1 ) 

YPLUS(I) 

UPLUS(I) 

HPLU$( I ) 

1 

O.OOOOE+OO 

O.OQOOE+QO 

0.13436+03 

0.00006+00 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

O.OOOOE+OO 

2 

0.2000E-02 

0.1111E+03 

0 . 1344E+03 

0.00006+00 

O.OOOOE+OO 

0.1272E+02 

0.5690E+01 

O.OOOOEtOO 

3 

0.6174E-02 

0.3774E+03 

0.13446*03 

0.00006+00 

O.OOOOE+OO 

0.3925E+02 

0.19336+02 

O.OOOOE+OO 

4 

0.6930E-02 

0 . 4100E+03 

0.13456+03 

0.00006+00 

O.OOOOE+OO 

0.4406E+02 

0.2100E+02 

O.OOOOE+OO 

5 

0.7686E-02 

0 . 4408E+03 

0.1344E+03 

O.OOOOE+OO 

O.OOOOE+OO 

0.48866+02 

0.2258E+02 

O.OOOOE+OO 

6 

0.8442E-02 

0.47046+03 

0.13446+03 

0.00006+00 

0.00006+00 

0.5367E+02 

0.2410E+02 

O.OOOOE+OO 

7 

0.9198E-02 

0 . 4989E+03 

0.13456+03 

0.00006+00 

O.OOOOE+OO 

0.58486+02 

0.25566+02 

O.OOOOE+OO 

8 

0.9955E-02 

0 . 5266E+03 

0.1344E+03 

0.00006+00 

O.OOOOE+OO 

0.63296+02 

0.26986+02 

O.OOOOE+OO 

9 

0.1071E-01 

0.5536E+03 

0.13446+03 

0.00006+00 

0.00006+00 

0.6809E+02 

0.2836E+02 

0.00006+00 

10 

0.1147E-01 

0.58016+03 

0.13446*03 

0.00006+00 

0.00006+00 

0.72906+02 

0.2972E+02 

O.OOOOE+OO 

11 

0 . 1222E-01 

0.6062E+03 

0.13446+03 

0.00006+00 

0.00006+00 

0.77716+02 

0.31056+02 

O.OOOOE+OO 

12 

0.1298E-01 

0.6319E+03 

0.13456+03 

O.OOOOE+OO 

O.OOOOE+OO 

0.8254E+02 

0.3237E+02 

O.OOOOE+OO 

13 

0.1373E-01 

0.65736+03 

0.13446+03 

0.00006+00 

0.00006+00 

0.87316+02 

0.3367E+02 

0.00006+00 

14 

0 . 1449E-01 

0 . 6826E+03 

0.13446+03 

0.00006+00 

0.00006+00 

0.9214E+02 

0.3497E+02 

O.OOOOE+OO 

15 

0.1524E-01 

0.70786+03 

0.13446+03 

O.OOOOE+OO 

0.00006+00 

0.96916+02 

0.3626E+02 

0.00006+00 

16 

0.1600E-01 

0.73306*03 

0.13446+03 

0,00006+00 

0.00006+00 

0.10176+03 

0.3755E+02 

O.OOOOE+OO 

17 

0.16766-01 

0.75826*03 

0.13456+03 

0.00006+00 

0.00006+00 

0.10666+03 

0.38846+02 

0.00006+00 

18 

0.1751E-01 

0.78356*03 

0.13456+03 

O.OOOOE+OO 

0.00006+00 

0.11136+03 

0.4014E+02 

0.00006+00 

19 

0.1827E-01 

0.80896*03 

0.13456+03 

0.00006+00 

0.00006+00 

0.11626+03 

0.41446+02 

0.00006+00 

20 

0.1903E-01 

0.83456*03 

0.13446+03 

O.OOOOE+OO 

0.00006+00 

0.1210E+03 

0.42756+02 

0.00006+00 

21 

0.1978E-01 

0.36046*03 

0.13456+03 

0.00006+00 

0.00006+00 

0.12586+03 

0.44086+02 

0.00006+00 

22 

0.2054E-01 

0.88676*03 

0.13446+03 

O.OOOOE+OO 

0.00006+00 

0.1306E+03 

0.4542E+02 

O.OOOOE+OO 

23 

0 . 2129E-01 

0.9133E*03 

0.13446+03 

O.OOOOE+OO 

0.00006+00 

0.13546+03 

0.46796+02 

0.00006+00 

24 

0.2205E-01 

0.94046*03 

0.13456+03 

0.00006+00 

0.00006+00 

0.14026+03 

0.48186+02 

0.00006+00 

25 

0.2281E-01 

0.96816+03 

0.13446+03 

0.00006+00 

0.00006+00 

0.14506+03 

0.49596+02 

O.OOOOE+OO 

26 

0.2356E-01 

0.99656+03 

0.13446+03 

0.00006+00 

O.OOOOE+OO 

0.1498E+03 

0.51056+02 

0.00006+00 

27 

0.2432E-01 

0.10266*04 

0.13446+03 

0.00006+00 

O.OOOOE+OO 

0.15466+03 

0.52566+02 

0.00006+00 

28 

0.2507E-01 

0.10566+04 

0.13456+03 

0.00006+00 

O.OOOOE+OO 

0.15946+03 

0.5410E+02 

O.OOOOE+OO 
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29  0.2583E-I1  0 .1C87E*C4  0.1  345E+03  0  .COCCE+OO 

30  0.2659E-01  0.:il9£+04  0.13446+03  0 .COOOE+OO 

31  0 .27346-01  0 .11526+04  0.1344E+03  O.OCOOE+OO 

32  0 . 28106-01  0 . 1 188E+04  0.1345E+03  0 .OOOOE+OO 

33  0.2886E-01  0.1225E+04  0.1345E+03  0. 0000E+00 

34  0.2961E-01  0 .1264E+04  0 .1344E+03  0. OOOOE+OO 

35  0.3037E-::  0 .13066+04  0.1344E+03  0. OOOOE+OO 

36  0 .31126-01  0 . 1 352E+04  0.1344E+03  0. OOOOE+OO 

37  3.3138E-01  0 .14026+04  0.1345E+03  0. OOOOE+OO 

38  C.3264E-01  0.1457E+04  0.1344E+03  0. OOOOE+OO 

39  3.3339E-01  0 . lr20E+04  0.1344E+03  O.OCOOE+OO 

40  0.33646-01  0 . 1544E+04  0.1344E+03  0. OOOOE+OO 

41  0.3390E-01  0 . 1 569E+04  0.1343E+03  C. OOOOE+OO 

42  3.34156-01  0.1597E+04  0.1342E+03  0. OOOOE+OO 

43  0.3440E-01  0.1626E+04  0.1341E+03  0. OOOOE+OO 

44  0.3465E-01  0.1657E+04  0.1339E+03  O.OCOOE+OO 

45  0.3491E-01  0.1688E+04  0.1337E+03  0, OOOOE+OO 

46  0 . 35166-01  0.1721E+04  0.1335E+03  0. OOOOE+OO 

47  C.3541E-01  0.1755E+04  0.1334E+03  0. OOOOE+OO 

48  C .35666-01  0.17896+04  0.13326+03  0. OOOOE+OO 

49  0.3592E-C1  0.1824E+04  0.1331E+03  0. OOOOE+OO 

50  0.3617E-01  0.1859E+04  0.1329E+03  0. OOOOE+OO 

51  0.3642E-01  0.1894E+04  0.1328E+03  0. OOOOE+OO 

52  0.3668E-01  0.1929E+04  0.13286+03  0. OOOOE+OO 

53  0 . 36936-01  0.1964E+04  0.1328E+03  0. OOOOE+OO 

54  0.3718E-01  0.1998E+04  0.1328E+03  0. OOOOE+OO 

55  0.3743E-01  0.2031E+04  0.1328E+03  0. OOOOE+OO 

56  0 .37696-01  0.2063E+04  0.1329E+03  0. OOOOE+OO 

57  0.3794E-01  0.2095E+04  0.1330E+03  O.OOOOE+OO 

58  0 . 38196-01  0.2124E+04  0.1331E+03  O.OOOOE+OO 

59  0.3844E-01  0.2153E+04  0.1333E+03  O.OOOOE+OO 

60  0.3870E-01  0.2179E+04  0.1334E+03  O.OOOOE+OO 

61  0 . 38956-01  0.22046+04  0.1336E+03  O.OOOOE+OO 

62  0.3920E-01  0.2226E+04  0.1337E+03  O.OOOOE+OO 

63  0.3945E-01  0.2246E+04  0.1338E+03  O.OOOOE+OO 

64  0 . 39716-01  0.2264E+04  0.1340E+03  O.OOOOE+OO 

65  0.3996E-01  0.2279E+04  0.1342E+03  O.OOOOE+OO 

66  0 . 40216-01  0.2290E+04  0.1342E+03  O.OOOOE+OO 

67  0 .40466-01  0.2299E+04  0.1343E+03  O.OOOOE+OO 

68  0.4072E-01  0.2305E+04  0.1344E+03  O.OOOOE+OO 

69  0 . 4097E-01  0.2307E+04  0.1343E+03  O.OOOOE+OO 

70  0.4097E-01  0.2306E+04  0.1344E+03  O.OOOOE+OO 


INTS  XU  USU  K  F  REM 

HE  EXX 

rHE  TURBULENT  COUETTE  WALL  FUNCTION  IS  BEIN6  USED,  AT  1 


O.OOOOE+OO  0.1642E+03  C.5S69E*02  0.3000E*00 
0 . 00006*00  0.1691E+03  0.5732E+02  O.OCOOE+CO 
O.OCOOE+OO  0.1738E+03  0.5902E+02  O.OOOOE+OO 
O.OOOOE+OO  0.1787E+03  0.6086E+02  O.OOOOE+OO 
O.OOOOE+OO  0.1835E+03  0.6275E+02  O.OOOOE+OO 
O.OOOOE+OO  0 . 1883E+03  0.6475E+02  O.OOOOE+OO 
O.OOOOE+OO  0.19316+03  0.669CE+02  O.OOOOE+OO 
O.OOOOE+OO  0 . 1979E+03  C.6926E+02  0  OOOOE+OO 
O.OCOOE+OO  0.2027E+03  0.71826+02  O.CCOOE*GO 
O.OOOOE+OO  0.20756*03  0.7464E+02  O.OOCCE*OC 
O.OCOOE+OO  0.2123E+03  0.7787E+02  O.OCCOE*CO 
O.OOOOE+OO  0.2139E+33  0.79106*02  O.OOOOE+OO 
O.OOOOE+OO  0.2155E+03  0.8038E+02  O.OOOOE*CC 
O.OOOOE+OO  0.2171E+03  0.8181E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2187E+03  0.8330E+02  O.OOOOE+CC 
O.OOOOE+OO  0.2203E+03  0.8489E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2220E+03  0.8647E+02  O.OOOOE+OO 
O.OOOOE+OO  0 . 2235E+03  0.8816E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2251E+03  0.8991E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2267E+03  0.9165E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2284E+03  0.9344E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2300E+03  0.9523E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2316E+03  0.9703E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2332E+03  0.9882E+02  O.OOOOE+OO 
O.OOOOE+OO  0.2348E+03  0.1006E+03  O.OOOOE+OO 
O.OOOOE+OO  0.23646+03  0.1024E+03  0.00006+00 
O.OOOOE+OO  0.2380E+03  0.1040E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2396E+03  0.1057E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2412E+03  0 .10736+03  O.OOOOE+OO 
O.OOOOE+OO  0.2428E+03  0.1088E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2444E+03  0.1103E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2460E+03  0.1116E+03  O.OOOOE+OO 
O.OOOOE+OO  0 .24766+03  0.11296+03  O.OOOOE+OO 
O.OOOOE+OO  0.2492E+03  0.1140E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2508E+03  0.1151E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2525E+03  0.U60E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2541E+03  0.1167E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2556E+03  0.11736+03  O.OOOOE+OO 
O.OOOOE+OO  0.2572E+03  0.1178E+03  O.OOOOE+OO 
O.OOOOE+OO  0.2589E+03  0.11816+03  0 .00006+00 
O.OOOOE+OO  0.2605E+03  0.11826+03  O.OOOOE+OO 
O.OOOOE+OO  0.2605E+03  C .11816+03  O.OOOOE+OO 


CF2  H  REH  ST  F(l,ilALL)  A 


EAST  AT  THIS  INTEGRATION 


ROUTINE  LANSU8  HAS  8EEN  CALLED 
N  HAS  SHIFTED  TO  68  AT  INTS  *  71 

=  81  0 .94056+00  0.2306E+04  O.OOOE+OO  O.OOOE+OO  0.5915E+04  0.10186-02  0.15116+01  0.54746+05  -.2737E-02  0.13436*03  -.34 
8E-01  -.138E+04 
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/  ’  ' 
''l .  / 


0. OOOOE+OO 

0. OOOOE+OO 

*) 

L. 

0.4689E-O3 

0.32072*03 

3 

0.2256E-02 

0.1136E*04 

4 

0 .38262*02 

0.124CE+04 

5 

0.45172*02 

0.12732*04 

6 

0.5246E-C2 

0 . 1 3C4E+04 

7 

0.6C08E-C2 

0.13322*04 

8 

0.;303E-02 

0.13592+04 

9 

0.76292*02 

0.13832*04 

•  A 
.  J 

0.34832*02 

0.14072+04 

.  4 

0 .93652*02 

0.14292*04 

12 

0.10272*01 

0.14512*04 

13 

a  +  *  a*  ; _a< 

V.iUti.  V* 

0.1472E-C4 

14 

0. 12162*01 

0.14932+04 

15 

0.131 5E-01 

0.15152+04 

16 

C.1415E-01 

0.15382+04 

17 

0.15182*01 

0.15622*04 

18 

0.16232*01 

0.1587E+04 

19 

0 . 1 729E-01 

0 . 16 12E+04 

20 

0.1837E-01 

0.1638E+04 

21 

0.1948E-01 

0.16642+04 

22 

0.2058E-01 

0.1691E+04 

23 

0.2172E-01 

0.1718E+04 

24 

0.2286E-C1 

0.1746E+04 

25 

0.2402E-01 

0.17742+04 

26 

0 .25212-01 

0.1803E+04 

27 

0.2639E-11 

0.1831E+04 

28 

0. 27612*01 

0.1859E+04 

29 

0. 28822*01 

0.18882+04 

30 

0.3007E-01 

0.19172+04 

31 

0.3134E-01 

0 . 1945E+04 

32 

0. 32612*01 

0.1974E+04 

33 

0.3391E-01 

0.2002E+04 

34 

0.3524E-01 

0.2030E+04 

35 

C.3657E-01 

0.2057E+04 

36 

0.3795E-01 

0.2084E+04 

37 

0.3934E-01 

0.211 1E+04 

38 

0. 4078E-01 

0.2137E+04 

39 

0.4225E-01 

0.2162E+04 

40 

0.4375E-01 

0.2186E+04 

41 

0.4426E-01 

0.2193E+04 

42 

0.4480E-01 

0.2201E+04 

43 

0.4532E-01 

0.2209E+04 

44 

0.4586E-01 

0.2216E+04 

45 

0.464CE-01 

0.2223E+04 

46 

0. 46972*01 

0.22312+04 

47 

0.4753E-01 

0.2238E+04 

48 

0. 48102*01 

0.22442+04 

49 

0.4868E-01 

0.2251E+04 

50 

0.4930E-01 

0.22572+04 

51 

0.4990E-01 

0.22632+04 

rr  •  • '  :'  *)  -  \ 

i  *  i  *  /  (  <  *  / 

0.1343E*C3  0.0C0<5E*00 
0.13432*03  0.00002*00 
0 . 1342E*03  0 .OOOOE+OO 
0 . I342E+03  0. OOOOE+OO 
0 . 1342E+03  0 .OOOOE+OO 
0 .13422+03  3. OOOOE+OO 
C.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0 .1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
C.1342E+C3  0. OOOOE+OO 
0.13422+03  0. OOOOE+OO 
0 . 1 242E+03  0. OOOOE+OO 
0.1342E+C3  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.13422+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+C3  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0 .OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0 . 1 342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1342E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0 . 1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0  OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0 . 1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0 . 1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0 . 1 343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 
0 . I 343E+03  0. OOOOE+OO 
0.1343E+03  0. OOOOE+OO 


:!3,:;  y°:jsc) 

0. OOOOE+OO  0. OOOOE+OO 
0. OOOOE+OO  0.1123E+02 
0. OOOOE+OO  0.5403E+02 
0. OOOOE+OO  0 . 9165E+02 
O.OCOCE-OO  0.10822*03 
0.00002*00  0.1257E+C3 
0.00002*00  0 . I439E*03 
0 ,OOOCE*CO  0 . 1 630E  +03 
0. OOOOE+OO  0.1827E+03 
0. OOOOE+OO  0.2032E+03 
0. OOOOE+OO  0.22432*03 
0. OOOOE+OO  0 .2461E+Q3 
0. OOOOE+OO  0.26862+03 
0. OOOOE+OO  0.2914E+03 
0. OOOOE+OO  0.3151E+03 
0. OOOOE+OO  0.3389E+03 
0. OOOOE+OO  0 . 3636E+03 
0. OOOOE+OO  0.3888E+03 
0. OOOOE+OO  0.4141E+03 
0. OOOOE+OO  0.4401E+03 
0. OOOOE+OO  0.4665E+Q3 
0. OOOOE+OO  0.4930E+03 
0. OOOOE+OO  0.5202E+03 
0. OOOOE+OO  0.5475E+03 
0. OOOOE+OO  0.5755E+03 
0. OOOOE+OO  0 . 6038E+03 
0. OOOOE+OO  Q  6322E+03 
0. OOOOE+OO  0.6613E+03 
0. OOOOE+OO  0.6904E+03 
0. OOOOE+OO  0.7203E+03 
0. OOOOE+OO  0.7507E+03 
0. OOOOE+OO  0.7811E+03 
0. OOOOE+OO  0.8124E+03 
0. OOOOE+OO  0.8442E+03 
0. OOOOE+OO  0.8761E+03 
0. OOOOE+OO  0.9091E+03 
0. OOOOE+OO  0.9423E+03 
0. OOOOE+OO  0.9768E+03 
0. OOOOE+OO  0 . 1012E+04 
0. OOOOE+OO  0.1048E+04 
0. OOOOE+OO  0 . 1060E+04 
0. OOOOE+OO  0.1073E+04 
0. OOOOE+OO  0.1086E+04 
0. OOOOE+OO  0.1098E+04 
0. OOOOE+OO  0.1111E+04 
0. OOOOE+OO  0.1125E+04 
0. OOOOE+OO  0.1139E+04 
0. OOOOE+OO  0.1152E+04 
0. OOOOE+OO  0.1166E+04 
0. OOOOE+OO  0.1181E+04 
0. OOOOE+OO  0.1195E+04 


0. OOOOE+OO  0.00002*00 
0.1116E*02  0.93512*01 
0.15452*02  0.13512*02 
0 . 1 686E+02  0.1473E+02 
0.17312+02  0.15032*02 
0 . 1 773E+02  0.15432*02 
0 . 181 1E+02  0.15912*02 
0.1847E*02  0 .15982*02 
0 . 1881E+02  0.16082*02 
0.1913E+02  0.16412*02 
0.1943E+02  0.16432*02 
0 ,1973E*C2  0 . 1683E*02 
0.2001E+C2  0. 16662+02 
0 . 203CE+02  0.16772*02 
0.2060E+02  0.1689E+02 
0.2091E-02  0.1709E+02 
0.2124E+02  0.1717E+02 
0.2157E+02  0.17242+02 
0.2191E+02  0.1712E+02 
0.2227E+02  0.1712E+02 
0.2263E+02  0.1692E+02 
0.2299E+C2  0.1687E+02 
0.233 7E+02  0.1662E+02 
0.2374E+02  0.1643E*02 
0.2412E+02  0.1599E+02 
0 .2451E+02  0.1559E+02 
0.2489E+02  0.1512E+02 
0.2528E+02  0.1453E*02 
0.2567E+02  0.1399E+02 
0.2606E+02  0.1330E+02 
0.2645E+02  0.1243E+02 
0.2683E+02  0.1160E+02 
0.2722E+02  0.1060E+C2 
0.2760E+02  0.9536E+01 
0.2797E+02  0.8306E+01 
0.2834E+02  0.7000E+01 
0.2870E+02  0.5587E+01 
0.2905E+02  0.4062E+01 
0.2939E+02  0.2467E+01 
0.2972E+02  0.8157E+0C 
0 . 2982E+02  0.1824E+00 
0.2993E+02  -.3913E+00 
0.3003E+02  -.9981E+00 
0.3013E+02  - .1645E+01 
0.3023E+02  - .2205E+01 
0.3033E+02  - .2881E+01 
0.3042E+02  - .  3438E  +  C 1 
0.3051E+02  - . 4085E *01 
0.3060E+02  -.4732E+01 
0.3069E+02  -.5375E+01 
0.3077E+02  -.60612*01 
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32 

0 . 305' E-91 

0.2269E*04 

53 

0.5116E-01 

0 . 2275E+04 

54 

0.5179E-01 

0.22806+04 

55 

0 .5243E-01 

0.2285E+04 

56 

0. 53086-01 

0.2289E+04 

0.53776-01 

0.2293E+04 

58 

0.5444E-01 

0.2297E+04 

59 

0 . 55116-01 

0.2299E+04 

60 

0.5580E-01 

0 . 23026+04 

61 

0.36326-01 

0.2304E+04 

62 

0.5722E-01 

3.2305E+C4 

63 

0 .57926-01 

0 .23066+04 

64 

0.5364E-01 

0 .23068+04 

65 

0 .5939E-01 

0 . 2306E+04 

66 

**  4A1  «  C./M 
v .Ov*«L  vi 

0 .23066+04 

67 

0.6083E-01 

0. 23068+04 

68 

0.6156E-01 

0.23C6E+04 

69 

0.6232E-01 

0 . 2306E+04 

70 

0.6305E-01 

0.23066+04 

71 

0.6305E-01 

0 . 2306E+04 

0.1343E*03  O.OCOCE*CO 
0.1343E*03  C . OOOCE+OO 
0.1343E-03  O.OOOOE+OO 
0.1343E*03  O.OOOOEtOO 
0.1343E*03  O.OOOOE+OO 
0 . I 343E+03  O.OOOOEtOO 
0 . 1343E+03  O.OOOOEtOO 
0.1343E+03  O.OOOOEtOO 
0 . 1 344E+03  O.OOOOEtOO 
0 .13446*03  O.OOOOEtOO 
0 . 1 344£t03  O.OOOOEtOO 
0.1344E+03  O.OOOOEtOO 
0 ,1344£t03  O.OOOOEtOO 
0.1344Et03  O.OOOOEtOO 
0 . 1 344£t03  O.OOOOEtOO 
0.1344E+03  O.OOOOEtOO 
0.1344E+03  0  OOE+OO 
0 .1344E+03  0 .UOUOE+OO 
0 . 1344£t03  O.OOOOEtOO 
0 . 1344£t03  O.OOOOEtOO 


O.OCOOE’OO  0.121CE*04 
O.OOOOE+OO  C.1225E+04 
O.OOOOE+OO  0.1241E+04 
O.OOOOEtOO  0.1256E+04 
O.OOOOEtOO  0 . 1271E+04 
O.OOOOE+OO  0 . 12886+04 
O.OOOOEtOO  0 . 1304E+04 
O.OOOOEtOO  0.I320E+04 
O.OCCOE+OO  0.13376+04 
0 . OOOOE+OO  0 . 1354E+04 
O.OOOOE+OO  0.1371E+04 
O.OOOOE+OO  0.13886+04 
O.OOOOEtOO  0.1405E+04 
O.OOOOE+OO  0.1423E+04 
O.OOOOE+OO  0 . 14406+04 
O.OOOOE+OO  0.1457E+04 
O.OOOOE+OO  0.1475E+04 
O.OOOOE+OO  0.1493E+04 
O.OOOOE+OO  0 . 1510E+04 
O.OOOOE+OO  0 . 1S10E+04 


0.3085E+02 

-,6668E*01 

0 . 3093E+C2 

- . 73216*01 

0.3100E+02 

- . 7931E+01 

0.3106E+C2 

-.85258+01 

0.31 12E+02 

- . 9069E+01 

0.31 18E+02 

- .96196+01 

0 . 31 22E+02 

-.1011E+02 

0 .3126E+02 

- . 1054E+02 

0.3130E+02 

-,1093E*02 

0 .31326+02 

- . 1 1 26E+02 

0.31348*02 

-.11486-02 

0.3135E+02 

-.11638*02 

0.3135E+02 

-.11668*02 

0.3135E+02 

- . 1 166E+02 

0.31356+02 

-.11658*02 

0 .31356+02 

- . 1 166E+C2 

0.3135E+02 

- . 1 166E+02 

0.3135E+02 

- .  1 166E+02 

0.3135E+02 

- . 1 166E+02 

0.3135E+02 

-.1166E+02 

The  oounflary  layer  thicknesses  at  xl( ft )=  0.9405 

Gelta99( in )=  0.6264  dispiaceient  deltal(in)-  0.3781  eoientui  delta2(in)*  0.0360 


74 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  NO.  0704-0188 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

AFIT /GA-'./iN  Y/S9  D-C4 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 

School  of  Engineering  AFIT/SNY 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 


7b.  ADDRESS  (C/ty,  State,  and  ZIP  Code) 


Air  Force  Institute  of  Technology 
■■.’right-Patterson  AF3f  Ohio  45433-6583 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Aright  Research  Develop  Ctr 


8c.  ADDRESS  (C/fy,  State,  and  ZIP  Code) 


8b.  OFFICE  SYM80L  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

?im 


10.  SOURCE  OF  FUNDING  NUMBERS 


'.•Jright-Patterson  AFB,  Ohio  45433 


PROGRAM  I  PROJECT 
ELEMENT  NO.  I  NO. 


TASK 

WORK  UNIT 

NO 

ACCESSION  NO. 

1 1 .  TITLE  (Include  Security  Classification) 

A  Computational  Model  For  Thickening  Boundary  Layers  With  Mass  Addition  For  Hypersonic 
engine  Inlet  Testing 


12.  PERSONAL  AUTHOR(S) 

Robert  D.  Clausen,  Japt,  USA? 


13a.  TYPE  OF  REPORT  13b.  TIME  COVERED  114.  DATE  OF  REPORT  (Year,  Month,  Day)  115.  PAGE  COUNT 

Thesis  from _ to _  I  1989  December  I  8“ 


16.  SUPPLEMENTARY  NOTATION 


17. _  COSATI  COOES _ I  18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

FIELD  I  group  |  sub-group  |  Boundary  Layers,  Inlets,  Mass  Addition,  Injection, 

Blowing 


19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Thesis  Advisor:  Lt.  Col.  Paul  I.  King,  Ph.D. 

Assistant  Professor  of  Aerospace  engineering 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

B  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  OTIC  USERS  UNCLASSIFIED 


DO  Form  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


A  computational  model  for  thickening  boundary  layers 
with  mass  addition  is  developed.  The  phenomena  of  uniform 
injection  into  a  two-dimensional  supersonic  stream  and 
subsequent  boundary  layer  growth  downstream  is  discussed. 
Analysis  of  the  injection  region  provides  the  thickness  of 
the  boundary  layer  just  aft  of  injection.  An  injection 
region  velocity  profile  is  then  used  to  approximate  the 
boundary  layer  profile  just  aft  of  injection  and  is  input 
into  a  finite-difference  boundary  layer  code.  Downstream 
profiles  and  thicknesses  are  calculated  and  compared  to 
experimental  results. 

The  computational  model  developed  here  provides  a  tool 
for  the  design  of  a  boundary  layer  generation  system  for 
hypersonic  engine  inlet  testing.  This  mass  addition  system 
is  needed  to  simulate  the  boundary  layer  developed  on  the 
forebody  of  hypersonic  vehicles.  An  example  is  discussed 
which  increases  the  natural  boundary  layer  thickness  17 
times . 


